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EFFECTS OF NON-LINEAR PARTITIONING 
BEHAVIOR ON NAPL CHARACTERIZATION 
VIA PARTITIONING TRACER TESTS 

by Elizabeth A. Fitzpatrick 

Chairperson of the Supervisory Committee: Professor William R. Wise 
Department of Environmental Engineering Sciences 

This report focuses on the non-linearities involved between tracer/non-aqueous phase 
liquid (NAPL) pairs and their subsequent effects on subsurface contaminant 
characterization. Two computer models were utilized; one model predicted the extent of 
the non-linear behavior between selected tracer/NAPL pairs. The second model utilized 
Freundlich and linear isotherm data fits to the non-linearities and predicted effluent 
response curves based on model conditions. Accepted methods for NAPL quantification 
in the subsurface were applied to the linear and non-linear effluent response curves and the 
differences were compared. 

The analyzed tracers include methanol, 2-methyl-2-hexanol, and 3-methyl-2-hexanol. 
Trichloroethylene (TCE) served as the NAPL for the computer simulations. 

The simulation results show that NAPL saturation calculations are systematically 
underestimated when linear tracer partitioning is assumed. The NAPL saturation was 
underestimated by as much as forty percent for the specific interactions and model 
conditions used during this investigation. 
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CHAPTER 1 

OVERVIEW 

Introduction 

Organic contaminants have become a problem in many locations throughout the country 
and the world. Over the past several decades, many sites have experienced leaks, spills, 
and the disposal of organic wastes. Many of the organic contaminants which exist at these 
sites are “non-aqueous phase liquids (NAPLs).” NAPL infiltration into the subsurface 
and contaminant transport via groundwater have caused vast expanses of hazardous 
waste sites. Efforts to remediate NAPLs from the subsurface are often hampered by the 
fact that it is difficult to characterize the amount and extent of the contaminant distribution 
in the ground. Conventional techniques have proven to be relatively ineffective for this 
purpose (Pope et. al., 1994). A recently developed technique to characterize NAPL 
distributions and amounts utilizes a technique which is typically used in oil field 
applications; this method involves the use of partitioning tracers. These tracers, when 
injected into the subsurface, partition into the NAPL phase; the amount of NAPL can be 
determined by measuring the tracer concentration downgradient from the tracer injection 
point. Several field demonstrations have already been completed to quantify residual 
NAPLs in the saturated zones of contaminated sites (Annable et al., 1994). Thus far, 
these demonstrations have assumed that a linear relationship exists between tracer/NAPL 
interactions. Specifically, the amount of tracer which partitions into the NAPL is 
determined by a single partition coefficient regardless of the relative amounts of tracer and 
NAPL present. This paper examines the combination of several already established 
computer modeling techniques to accomplish the following: (1) determine the extent of 
tracer/NAPL non-linearities, and (2) estimate the non-linear effects on NAPL volume 
calculations using simulation techniques. The modeling programs used during this study 
include the UNEFAC (Chen, 1992) and COMB ALL (Wise, 1985) models. 

Background 

Current methods used to characterize subsurface NAPL concentrations include well 
installation, core sampling and analysis, and borehole geophysical techniques (Pope etal., 
1994). These technologies are often costly, time consuming, and cause soil disturbance. 
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In addition, they cannot accurately predict NAPL distributions because samples are 
evaluated at distinct locations throughout the site. The presence of both soil 
heterogeneities and anisotropic conditions often lead to inaccurate NAPL characterization. 
As a result, remediation technologies employed at many sites cannot effectively remove 
the targeted contaminants. Thus, costs and remediation times are often underestimated. 

A recently developed technique to characterize NAPL concentrations and distributions 
employs the use of partitioning tracers. Partitioning tracers have been used successfully in 
oil field applications (Pope et al., 1994). In NAPL characterization, laboratory 
experiments and field demonstrations have already been conducted to determine 
partitioning tracer effectiveness (Annable et al., 1994; Pope et al., 1994). These 
applications are based on the knowledge that injected tracers partition into and out of the 
oil (NAPL) phase as they move through an aquifer; consequently, these tracers experience 
retardation. If the relationship between each tracer and the NAPL is known, then the 
amount of NAPL present along a streamline in an aquifer can be estimated based on the 
retardation of each tracer. Previous studies to determine NAPL presence in saturated 
media have included the use of a “non-partitioning” tracer and several alcohols for 
partitioning tracers (Annable et al., 1994; Pope et al., 1994). The difference between the 
first moment of the non-partitioning and the partitioning tracer effluent curves can then be 
used to quantify certain NAPLS in the aquifer. 

Several advantages exist with the use of tracers. First, they only disturb the soil at tracer 
injection and extraction points. Second, they account for soil heterogeneities and 
anisotropic conditions by moving with the groundwater through the aquifer. Finally, only 
small amounts of tracers are required to adequately characterize the NAPL profiles. 
However, several limitations exist. One major limitation is that non-linearities exist 
between tracer/NAPL interactions. The amount of non-linearity depends on the mole 
fraction ratio between the tracer/NAPL pair. Therefore, the amount of retardation 
experienced by the partitioning tracer varies with the tracer/NAPL composition. 
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The laboratory and field experiments thus far have been analyzed under the assumption 
that the tracer/NAPL interactions are linear. However, additional information is required 
to fully understand the effects of non-linearities on tracer/NAPL interactions. The initial 
use of computer modeling techniques can predict the magnitude of non-linearities between 
tracer/NAPL pairs. This modeling may ultimately allow for more accurate contaminant 
level predictions by understanding the deviations from ideal conditions. There are, 
however, several disadvantages which exist to computer modeling techniques. First, the 
tracer/NAPL non-linearities are predicted by a computer model. Thus, these predictions 
may not be completely accurate. Second, the type of NAPL contaminant is not always 
known prior to characterization. The established computer modeling techniques can only 
predict non-linearities if the NAPL type is already known. 

Purpose 

This paper examines the extent of non-linearities involved between tracer/NAPL pairs and 
their effect on NAPL characterization. To accomplish this objective, tracer/NAPL non- 
linearities were predicted using the UNIFAC model (Chen, 1992). Established isotherm 
relationships were then “fit” to the UNIFAC predicted tracer/NAPL relationships. The 
isotherm fits, combined with hydraulic data, were used to determine partitioning tracer 
effluent curves through the use of the COMB ALL model (Wise, 1985). Finally, the 
COMB ALL effluent curve data was utilized to compare accepted linear methods for 
NAPL quantification to the actual non-linear behavior which actually exists. 
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CHAPTER 2 

MODELING PROCEDURE 
Model Parameters 

Three tracers were analyzed with one NAPL during this investigation. The analyzed 
tracers include methanol, 2-methyl-2-hexanol, and 3-methyl-2-hexanol. Methanol has 
often been used as a “non-reactive” tracer, while the other two represent typical 
partitioning tracers (Pope et al ., 1994). These tracers interacted with the NAPL 
trichloroethylene (TCE) for the purposes of this study. Table I provides the structural 
chemical makeup for each tracer and TCE.. 


Table I. Structural Formulas of Analyzed Compounds 



Trichloroethylene 

Methanol 

2- Methyl-2-Hexanol 

3- Methyl-2-Hexanol 


CIHC=Ca 2 

CH 3 OH 

CH 3 -CCH 30 H-CH 2 -CH 2 -CHrCH 3 

CH3-CHCH3-CHOH-CH2CH2-CH3 


The simulations were only conducted between binary systems; i.e., the interactions 
between one tracer and TCE were observed for each model simulation. In addition, the 
equations used in the model simulations were conducted with several dimensionless 
variables. These variables include the following: (1) Concentration, c; (2) Time, t; (3) 
Velocity, v; (4) Distance, xl, and (5) the Dispersion Coefficient, Dl. The models simulate 
a theoretical soil column in which no gas phase exists. Table II summarizes the 
parameters which were used throughout all model simulations. 
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Table II. Parameters Used for UNIFAC and COMBALL Simulations 


PORE 

POROSITY 

DISPERSION 

TEMP 

WATER 

NAPL 

INJECTION 

VOLUME 


COEFF. 


SAT. 

SAT 

TIME 

(PV) 

Ol) 

(DO 

(T) 

(Sw) 

(Sn) 

(At)* 

1.0 

30% 

0.01 

298K 

95% 

5% 

ilk:--'- • • • V" : ilhik • = " 

0.5 


*Injection time is expressed in terms of pore volumes for model simulations. 


Assumptions 

The UNIFAC and COMBALL runs were conducted under several general assumptions. 
First, the system is assumed to be at equilibrium. Second, the tracer only interacts with 
the NAPL in the column; i.e., no other tracer reactions such as biological or radioactive 
decay occur within the column. Third, no water phase exists within the tracer/NAPL 
body. Finally, the NAPL component is sparingly soluble in H 2 O. In other words, the 
entire NAPL volume in “the column” is available for tracer partitioning and accurate 
NAPL quantification. 


The UNIFAC Model 

The UNIFAC model predicts the interactions between each tracer/NAPL pair. Formally 
termed “The UNIFAC Group-Contribution Method,” this model was originally developed 
for chemical engineering purposes (Chen, 1992). The “group contribution” method 
assumes that compounds may be broken down into functional groups which have a unique 
contribution toward the compound property. By using this method, the UNIFAC model 
can predict liquid-phase activity coefficients. These activity coefficients represent the non¬ 
ideal behavior of a solution. The activity coefficient calculations for each molecule, /', are 
based on the following equation: 


(1) lny, =lny C j + lny 1 *\ 
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where Yi is the activity coefficient for the component, /, yi C is the combinatorial part, and 
Yi R is the is the residual part. The combinatorial portion takes into account the entropy 
effects and depends on the size and the shape of the molecule /. In the combinatorial part, 
the group surface volume and area of each functional group are used as model 
parameters. Chen (1992) presents the equations used to determine the ultimate activity 
coefficients for mixtures. 

The UNIFAC model was specifically utilized during this investigation to determine the 
activities between each tracer/TCE pair at varying tracer mole fractions between 0 and 
1.0. The proper use of UNIFAC requires user identification of the functional subgroups 
for each compound. The chemical configuration of each compound determines the 
functional subgroup types and subgroup quantities for that compound. The subgroup 
information for each compound must be entered at a specified reaction temperature and 
mole fraction. The UNIFAC program, with this information, calculates activity 
coefficients for each tracer/NAPL interaction between mole fractions of 0 to 1.0. Chen 
(1992) provides a table of functional subgroups and their “assigned” UNIFAC subgroup 
numbers. Table III provides the functional subgroup information for each tracer/TCE 


pair. 
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Table HI. UNIFAC Functional Subgroup Information for Tracers and NAPL 


[COMPOUND 

FUNC. 

MAIN 

SUB- 

TOTAL 


GROUP 

GROUP # 

GROUP# 

GROUPS 

Trichloroethylene 

H-C=C 

2 

””"Z’ 

. 1. 

E 

CI-(C=C) 

37 

70 

3 

Methanol 

CHjOH 

6 

16 

1 isi 

2-Methyl-2-HexanoI 

ch 3 

1 

1 

3 

p'v ; . 

ch 2 

1 

:- 2 :!: : 

3 


c 

1 

4 

1 

l:;; ; ; ; ; •’ :=. •. •’ 

OH 


15 


| 3-Methyl-2-Hexanol 

ch 3 

1 

1 

3 


ch 2 

. 1 

2 : 

2 

i 

CH 

1 

3 

2 

i-il" ■ ■■ t 

i 

OH 

; . 5 : . 

15 

1 


Appendix (A) includes an example UNIFAC input and output data file for the 2-methyl-2- 
hexanol/TCE system. The input file varies tracer mole fractions from 0 to 1.0 with an 
increasing mole fraction interval of 0.1. The corresponding TCE mole fractions range 
from 1.0 to 0 with a decreasing mole fraction interval of 0.1. Additional input parameters 
are further explained in Appendix (A). 

The UNIFAC output data files were then used to calculate the tracer and TCE activities; 
each activity coefficient multiplied by it’s corresponding mole fraction yields that mole 
fraction’s activity. Figures (1), (2), and (3) display the results for methanol/TCE, 2- 
methyl-2-hexanol/TCE, and 3-methyl-2-hexanol/TCE combinations. The lower “X” axis 
represents the tracer mole fraction, while the upper “X” axis represents the TCE mole 
fraction. 
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METHANOL/TRICHLOROETHYLENE SYSTEM 


MOLE FRACTION TCE 

1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 



MOLE FRACTION METHANOL 


-£-METHANOL ACTIVITY 

(ai) 

-METHANOL RAOULT'S 

LAW 

-D-TCE ACTIVITY (ai) 

-TCE RAOULT’S LAW 


Figure 1. 


2-METHYL-2-HEXANOL/TRICHLOROETHYLENE SYSTEM 


MOLE FRACTION TCE 

1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 



-£-2-METHYL-2-HEXANOL 

ACTIVITY (ai) 

-2-METHYL-2-H EXANOL 

RAOULTS LAW 
-□-TCE ACTIVITY (ai) 

-TCE RAOULTS LAW 


Figure 2, 
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3-METHYL-2-HEXANOL/TRICHLOROETHYLENE SYSTEM 


MOLE FRACTION TCE 

1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 



—6— 3-METHYL-2-HEXANOL 
ACTIVITY (ai) 

-3-METHYL-2-HEXANOL 

RAOULTS LAW 
—□— TCE ACTIVITY (ai) 


-TCE RAOULTS LAW 


Figure 3. 


The UNIFAC simulations represented by figures (1), (2), and (3) were conducted under 
the assumption that Henry’s Law is valid for both the tracer constituent and TCE; 
therefore, the activity in each figure represents the aqueous phase concentration of the 
tracer or TCE. Here, the aqueous phase concentration is defined as the compound’s 
actual concentration divided by the compound’s solubility. Figure (1), which represents 
the methanol/TCE interactions, displays highly “non-linear” results. Since methanol is 
highly soluble in water, it may not follow the Henry’s Law assumption. Thus, the 
methanol/TCE system requires further investigation to accurately quantify the portion of 
the non-linearities which result from TCE interactions. 


Figures (1), (2), and (3) also depict Raoult’s Law for these tracer/TCE systems; the 
aqueous phase concentration of the specified compound is given by its mole fraction since 
Yi (the activity coefficient) equals 1.0. The “non-ideal” behavior of each tracer/NAPL pair 
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can be easily seen by examining the UNIFAC results against Raoult’s Law. In all three 
cases, the non-idealities result in higher aqueous phase tracer concentrations than those 
predicted by Raoult’s Law. As mandated by thermodynamics, the compounds obey 
Raoult’s Law at mole fractions close to 1.0. 

Recent experiments have used relatively small concentrations of tracer to characterize 
NAPL concentrations (Annable et al., 1994; Pope etal., 1994). Therefore, the lower left 
portion of each of the preceding figures represents the primary region of interest; the non- 
linearities which exist in the lower left region must be further analyzed. For this reason, 
more detailed UNIFAC runs were conducted on smaller tracer mole fraction intervals; 
these tracer mole fraction intervals first ranged from 0 to 0.2 and then from 0 to 0.05. The 
activity (aqueous phase concentration) vs. tracer mole fraction in each case was plotted 
for analysis in figures (4) through (9). 


*Note: The notation “c” used in the text is represented by “Cw” in the UNIFAC figures. 




TRACER ACTIVITY VS. MOLE FRACTION FOR METHANOL/TCE 



TRACER ACTIVITY VS. MOLE FRACTION FOR METHANOL/TCE 

(WHEN USED AS A TRACER FOR TCE NAPL) 
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TRACER ACTIVITY VS. MOLE FRACTION FOR 2-METHYL-2-HEXANOL/TCE 



Figure 6. 


UNIFAC ACTIVITY VS. MOLE FRACTION FOR 2-METHYL-2-HEXANOL/TCE 



Figure 7. 















Figures (4), (6), and (8) demonstrate that tracer mole fraction intervals between 0 and 0.2 
display significantly non-linear behavior. The non-linear behavior is also evident between 
mole fraction intervals of 0 and 0.05. However, the non-linearities are not as pronounced 
within this smaller interval. 


Freundlich Relationship 

Each of the six “Aqueous Phase Concentration vs. Tracer Mole Fraction” plots was 
modeled to fit a Freundlich isotherm type relationship. For each tracer/TCE pair, the mole 
fraction interval which demonstrated a closer fit to the Freundlich model was chosen for 
further investigation. The following expression was approximated using data points from 
the UNIF AC-generated c-x relationship: 


(2) x=Pc“ 

x represents the tracer mole fraction 


moles tracer 


moles_tracer + moles_ TCE 

c represents the aqueous phase concentration (mass in solution/solution volume) 
[ML* 3 ] 

a represents the Freundlich exponent 
P represents the Freundlich coefficient 


The Freundlich parameters a and P were determined by minimizing the error between the 
UNIF AC calculated “x-c” relationship and the Freundlich approximated “x-c” relationship; 
a spreadsheet solver routine (Microsoft® Excel Solver) was used to vary a and P until the 
error between the data and the Freundlich approximation was minimized. Figures (10) 
through (15) display the Freundlich fits for each tracer/TCE pair for mole fraction 
intervals between 0 to 0.2 and 0 to 0.05. These figures plot the dependent variable, c, as 
the abscissa so that the shape of the approximated Freundlich relation can be easily 
recognized. 



TRACER MOLE FRACTlOh TRACER MOLE FRACTlOh 


15 


METHANOL/TCE UNIFAC RESULTS COMPARED WITH FREUNDLICH APPROXIMATION 



■ MOLE FRACTION TRACER 
CALC. FREUND. MOLE FRACTION 


FOR MOLE FRACTIONS 
BETWEEN 

0 —> 0.2 
BETA= 0.29323377 
ALPHA=2.50636922 


X = BETA*Cw A ALPHA 


Figure 10. 


METHANOL/TCE UNIFAC RESULTS COMPARED WITH FREUNDLICH APPROXIMATION 



-MOLE FRACTION TRACER 

.CALC. FREUND MOLE FRACTION 


FOR MOLE FRACTIONS 
BETWEEN 

0—> 0.05 
BETA= 0.13059633 
ALPHA=1.25733495 


X = BETA*Cw A ALPHA 

Figure 11. 









TRACER MOLE FRACTlOh TRACER MOLE FRACTK» 


2-METHYL-2-HEXANOL/TCE UNIFAC RESULTS COMPARED WITH FREUNDUCH 

APPROXIMATION 



-MOLE FRACTION TRACER 

- • - CALC. FREUND. MOLE FRACTION 

FOR MOLE FRACTIONS 
BETWEEN 

0 —> 0.2 
BETA= 1.26225144 
ALPA=1.92063211 


X ® BETA*Cw A ALPHA 


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Cw 


Figure 12. 


2-METHYL-2-HEXANOL/TCE UNIFAC RESULTS COMPARED WITH FREUNDUCH 

APPROXIMATION 


0.05 

0.045 

0.04 

0035 

0.03 

0.025 

0.02 

0.015 

0.01 

0.005 

0 



0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 


-MOLE FRACTION TRACER 

■ - • - CALC. FREUND. MOLE FRACTION 

FOR MOLE FRACTIONS 
BETWEEN 

0—> 0.05 
BETA= 0.41835494 
ALPHA=1.24654054 


X = BETA*Cw A ALPHA 


Figure 13. 
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3-METHYL-2-HEXANOL/TCE UNIFAC RESULTS COMPARED WITH FREUNDUCH 

APPROXIMATION 



-MOLE FRACTION TRACER 

.CALC. FREUND. MOLE FRACTION 


FOR MOLE FRACTIONS 
BETWEEN 

0 —> 0.2 
BETA= 1.24395432 
ALPHA=1.91242923 


X = BETA*Cw*ALPHA 


Figure 14. 


3-METHYL-2-HEXANOL/TCE UNIFAC RESULTS COMPARED WITH FREUNDUCH 

APPROXIMATION 



-MOLE FRACTION TRACER 
•CALC. FREUND. MOLE FRACTION 


FOR MOLE FRACTIONS 
BETWEEN 

0—> 0.05 
BETA= 0.41668518 
ALPHA=1.24437431 


X = BETA*Cw*ALPHA 


0 


Cw 


0.12 0.14 0.16 0.18 


Figure 15, 



















































































































































































































































































18 


In all three tracer cases, the Freundlich isotherms fit the UNIFAC data better at lower 
mole fraction intervals (x=0 to 0.05). The root mean-squared (RMS) approximation 
errors ranged from 5.112 X 10" 4 to 5.683 X 10' 4 for the tracers analyzed within the x=0 to 
0.05 interval (see Table IV). The higher mole fraction intervals (x=0 to 0.20) did not 
display a UNIFAC/Freundlich correlation as favorable as that for the lower mole fraction 
intervals. The RMS error ranged from 4.330 X 10' 3 to 8.846 X 10‘ 3 for the tracers 
analyzed within the x= 0 to 0.20 interval. Generally, the Freundlich calculated mole 
fractions fell below the UNIFAC fits at the lower portion of the mole fraction interval 
analyzed, while the Freundlich approximation resided above the UNIFAC results for the 
upper end of the analyzed mole fraction interval. The Freundlich approximation for the 
mole fraction interval x=0 to 0.05 was chosen for further analysis because of its close 
approximation to the UNIFAC results within this interval. Appendix (B) includes a 
sample spreadsheet solution for 2-methyl-2-hexanol/TCE (for x=0 to 0.05) and the 
method used to calculate the RMS error. 


Table IV. Tracer “3” and “a” Values With Error Approximations 
With Associated RMS Errors 


Tracer/TCE Pair 

x Interval 

3 

a 

RMS Error 

Methanol 

. 0^ 0.2 

.29324 

2.50637 

8.8458 X 10‘ 3 


0-» 0.05 

.13060 

1.25733 

5.6830 X 10* 4 

2-Methyl-2-Hexanol 

0^ 0.2 

1.26225 

1.92063 

4.3418 X IO" 3 

| 

j 

0-» 0.05 

.41836 

1.24654 

5.1531 X 10* 4 

i i^MethyW^iiandl 

0-> 0.2 

1.24395 

1.91243 

4.3298 X 10' 3 

1 

! 

L.-....... 

0^ 0.05 

.41669 

1.24437 

5.1122 X 10" 4 


Relation Between Mole Fraction (x) and Sorbed Concentration (q) 

Although the equation x=pc“ represents a Freundlich type equation, it must be related to 

the “classic” Freundlich isotherm so that its meaning can be fully understood. The 

“classic” Freundlich isotherm may be expressed in terms of sorbed concentrations: 
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(3) q=Bc A 

q represents the sorbed concentration (mass sorbed/solution volume) [ML' 3 ], 
c represents the aqueous phase concentration (mass in solution/solution volume) 
[ML' 3 ]. 

“B” and “A” represent Freundlich parameters. 

The relationship between the “classic” Freundlich isotherm and that which was solved 
during this investigation can be expressed using the following general analysis: 

(4) x= pc“ 

as determined by the approximations to the UNIFAC results (same as equation (2), 
but repeated here for complete analysis). Assuming that only one tracer exists in 
the “column,” 

, v moles tracer 

(5) x~— - - 

moles tracer + moles_NAPL' 

Further, assuming no gas phase is present, 

(6) Sn'Hiw - 1 

where s„ and Sw represent the NAPL (oil) and aqueous phase saturations, 
respectively. The tracer mass may be related to the water volume present with the 
following equation: 

_ mass tracer in oil_per bulk volume 
q s w 

where q [ML' 3 ] represents the local storage of the tracer in the oil (NAPL) phase. 

Assuming that the tracer’s presence does not significantly affect the water 
saturation, equation (7) may be written using equation (5): 
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where M.W.t represents the molecular weight of the tracer [M/mole]. 


The tracer’s presence should not significantly affect the water saturation, Sw, since 
Sw is significantly higher than s„(95% compared to 5%). Thus, even a significant 
change in Sn (i.e., 20%) will only result in a minute change in Sw (1.05%). 

mn„ 

(9) moles NAPL per bulk volume = — -. 

- M.W. n 


tj represents porosity 


void volume 
bulk volume 


p n represents the NAPL density [ML' 3 ]. 

M.W.„ represents the molecular weight of the NAPL [M/mole]. 


Letting 

(10) £ = — —- {moles _ NAPL per bulk _ volume ), 

S w 


Combining equations (4) and (11), 

02 ) 

The solution represented by equation (12) does not yield a “simple” Freundlich 
relationship which correlates “B” and “P” or “A” and “a.” Therefore, the 
spreadsheet solver routine was used to solve for B and A based upon equations (3) 
and (11). 

( 13 ) 
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Letting 

B 

(14) B’=—, 

X 

(15) B’c a= 7I^ 

Letting v|/=B’c A , the Microsoft® Excel spreadsheet solver routine was used to calculate a 
>” value for each x. Specifically, the routine varied the values of B’ and A until the 

x 

minimum error was reached between vi/ and -. Finally, “B” was found by evaluating 

1-x 

B’ with appropriate values of E, for each tracer/TCE pair. The TCE molecular weight and 
density used for the E, calculations were 131.5 g/mole and 1.46 g/ml, respectively. Table 
V provides the values used to determine the Freundlich parameters for each tracer/TCE 
pair. 


Table V. Values Used to Determine Freundlich Isotherm Parameters 


Tracer ? 

11.W.T S, (g/ml) B- 

r/mole) 

: . •. f ; ;V: ; ' 

B 

A : RMS 

Error 

Methanol 

32 .0058622 .14143 

.82911 

1.29760 6.67 X I O' 4 

2-Methyl-2-HexanoI 

116 .0212505 .47007 

9.98927 

1.28623 6.14 X 10” 4 

3-Methyl-2-HexanoI 

116 .0212505 .46806 

9.94644 

1.28392 6.10 X10" 4 


Appendix (C) contains a sample spreadsheet solution for 2-methyl-2-hexanol/TCE which 
solves for B and A based on P, a, and E,. 
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CHAPTER 3 
COMB ALL MODEL 

Overview 

The COMBALL program, which is written in FORTRAN, determines the effluent 
response curve for general multi-well injection tests (Wise, 1985). For the purposes of 
this investigation, the model solved the streamline advective transport equation for a soil 
column into which a pulse input of a sorbing tracer was injected at one end. Although 
COMBALL can determine the effluent response curves for sorption which follows either 
the Langmuir or Freundlich isotherm, the routine only modeled the effluent response curve 
for the unfavorable Freundlich isotherm during this investigation. COMBALL solves the 
streamline advective transport equation which corresponds to the Freundlich isotherm 
using the method of characteristics. This solution, coupled with an appropriate 
description of the system’s hydraulics, is utilized by the COMBALL model to determine 
the effluent concentration through time. The COMBALL program is attached as 
Appendix (D). 


Advective Transport Equation 

The advective transport equation is critical in the determination of the effluent response 
curve for a given tracer/NAPL pair. For an arbitrary partitioning tracer, the general form 
of the transport equation is represented by the following (after Wise, 1985; Rainwater et 
al ., 1987; Wise and Charbeneau, 1994): 

.... 3} dc & 

(16) a + a + * = 0 


dq/dt represents the temporal change in tracer storage in the NAPL (immobile) 
phase [ML- 3 r 1 ]. 

dddt represents the temporal change in tracer storage in water [ML* 3 !” 1 ]. 
dc/dx represents the divergence of the tracer advective flux through the colum. t 
represents travel time [T] and expresses distance in the form of time 
(t=xl/Vp); it relates the pore water velocity, v p [LT 1 ], and the time it takes 
the tracer to travel along a streamline to a distance “xl” along the 
streamline [L], x is used so that the hydraulics and chemistry of the 
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advective transport equation can be treated in an efficient and consistent 
manner. 


(17) 


COMBALL is predicated on the local equilibrium assumption (LEA) and on 
steady hydraulics throughout the test. Based on the LEA, 


dt 

and 


(#) 

( dc' 

\cb) 

l ~du 


(18) 


dj'] dc dc 

1+ *Ja + *"° 


This version of the transport equation demonstrates its quasi-linear, hyperbolic nature. 

Characteristic Solution to the Advective Transport Equation 
Utilizing the Method of Characteristics, the following base characteristics apply at the 

extraction well, and can be represented mathematically in the x-t plane (Wise, 1985): 

dt dq . 

(19) — = 1 + — isochore. 

dr c dc 

The “isochore” characteristic equation applies for paths of constant concentration leaving 
the t and x axes. Since characteristics are not allowed to cross, a “shock” condition 
results whenever characteristics meet. These concentration “jumps” in the x-t plane satisfy 
the following: 


dt , 

( 2 °) *1 


STEP 


A q 

= 1+ A 

Ac 


shock. 


Using the Freundlich relation in equation (3), equations (19) and (20) become: 
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( 21 ) 


dt 


dr 


\ + ABc ( a ~ x) 


isochore. 


(22) — =\ + Bc (Al) shock 

STEP 

The shock equation represents the curve in the x-t plane in which the wave steps from a 
composition of (c, Bc A ) to (0, 0). This condition will be met for all effluent results 
presented in this investigation. 


Application of Characteristic Solution to Streamline Transport Problem 
The three tracers which were analyzed in this investigation were fitted to unfavorable 

Freundlich isotherms. Therefore, this investigation made use of the method of 

characteristics solution through the COMBALL subroutine for an unfavorable isotherm 

(see Appendix (D)). For an injection time of At, the shape of the characteristic curve for 

an unfavorable isotherm (A>1) resembles figure (16): 
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For A>1, seven equations govern the characteristic solution. A description of the 
equations and the areas to which they apply are the following: 

As can be seen in figure (16), the endpoint (x p i a t, t p i at ) represents the intersection of the 
shock wave of the Co and c=Co isochores. These points are described by equations (23) 
and (24): 


At 

( )Tpl *-Bc 0 {A -' ] (A- 1) 

(\ +ABc n (A ~ x) )At 
(24) ‘" le, 


Between the c=0 and the c=Co isochores, the following equation provides the 
concentration at any point (x, t): 


(25) 




The following equations determine the x values at which an isochrone (line of 
constant time) passes the c=Co and the c=0 isochores, respectively: 


(26) x 0 (t)=t 


(27) r Co (t) = 


t 

(l + ABc 0 (A ~ x) ) 


XsiEp(t) is defined as the x value at which an isochrone, t, crosses the step 
concentration and can be found by the following: 
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(28) For t<t p i at : XsiEp(t)- (J^y 

1 +DC Q 

(29) For t>t p iat: XsTEp(t) results from solving equation (29.a) using Newton’s Method. 
Equation (29.b) results. 


(29.a) G(x) 


(29.b) G’(x)= 


(±i) fill r±di 

*=t { a) -tr [Aj +I [aJ 




(29. c) 




04 - 1 ) 


Again, the preceding equations are those which are utilized in the COMB ALL model for 
an unfavorable Freundlich isotherm. 


Fractional Breakthrough Curve 

A “theoretical” fractional breakthrough curve was used to determine the effluent results 
from the COMBALL routine. The solution to the advection-dispersion (A-D) equation 
was used to generate the fractional breakthrough curve for the simulated, one¬ 
dimensional, laboratory column experiment. The A-D equation for a pulse input of time 
At is the following: 


1 L-vt v r L L + vJ 

(30) F = ~ [eiffl —j— ) + exp(— )erfc(- 


"2y[Dj 


Dr 


2 jDj 


(Fetter, p.458). 


F represents the decimal fraction of the injected tracer concentration observed at 
the end of the column. 

L represents the column length [L], 

v x represents the average linear pore water velocity [LT 1 ]. 

t represents the time after tracer injection [T], 
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Dl represents the dispersion coefficient [L 2 T']. 

Since the second “erfc” term quickly approaches zero through time, the A-D equation can 
be approximated as the following: 


(31) F = - [erfc( 




)] 


Equation (31) describes the system hydraulics for a column of length “L.” Any distance 
which is less than L can now be expressed in terms of pore volume; the length of the 
column equals 1 pore volume. Remembering that xl = v x t, the term x replaces the “v x t” 
term by taking the velocity through the column as one pore volume per pore volume 
injected. Thus, equation (31) can be expressed as the following: 


(32) V=-erfc\ 


1 — T 






For a dispersion coefficient of .01 


pore _ volumes 2 
pore volume 


the fractional breakthrough curve 


shown in figure (17) results. 
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A-D SOLUTION FOR FRACTIONAL BREAKTHROUGH CURVE 
PORE VOLUMES, D=.01 



COMB ALL Simulations 

Simulations were conducted for methanol, 2-methyl-2-hexanol, and 3-methyl-2-hexanol. 
TCE was again used as the NAPL. The Freundlich coefficients based on the UNIFAC 
results were used as model inputs. As was stated in the “UNIFAC” portion of this paper, 
the tracer mole fraction interval between 0 and 0.05 was used for each simulation since 
there was minimal error between the UNIFAC results and the Freundlich solution. 
Simulations were completed with a maximum time of 10 and a time interval of .025 
between data points. In addition, all COMBALL simulations utilized the same fractional 
breakthrough curve shown in figure (17). Linear approximations and sensitivity analyses 
were then conducted for the tracers. A sample COMBALL input and output file for 2- 
methyl-2-hexanol/TCE is attached as Appendix (E). 

The initial runs were conducted using Co values corresponding to mole fractions of 0.025 
(as determined by UNIFAC). Figure (18) displays the results for the three tracers plus a 
theoretical “non-reactive” tracer. Table VI provides a summary of the parameters which 
were entered into the COMBALL data file for a tracer mole fraction of 0.025. 
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TRACER EFFLUENT CURVES USING FREUNDLICH APPROXIMATION 

FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 


ft 



Figure 18. 


Table VI. Freundlich Parameters for Tracer/NAPL Mole Fraction = 0.025 


Tracer/TCE Pair B 

A 

:■ I:::::;:;:::::::::;:::;: A 


^10H M X\€&ctivc v 

. 1 : : : .. 1 

Methanol 82911 

1.29760 .27262 j 

2-Methyl-2-HexanoI 9.98927 

1.28623 .10581 

3-Methyl-2-Hexanol 9.94644 

1.28392 .10572 


The tracer effluent curves demonstrate that 2-methyl-2-hexanol and 3-methyl-2-hexanol 
are significantly retarded by the TCE as compared to the methanol. Since the only 
structural difference between 2-methyl-2-hexanol and 3-methyl-2-hexanol is the location 
of the methyl group, they perform almost identically in the presence of TCE. The fact that 
their curves are indistinguishable in Figure (18) verifies this fact. As was stated in the 
“Procedure” portion of this paper, methanol is often used as the “non-reactive” tracer. 
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Compounds which are highly polar are often used as non-reactive tracers since they do not 
tend to partition into the NAPL. Non-polar tracers, however, partition into the NAPL 
more easily and have higher retardation factors. Figure (18) demonstrates that although 
methanol is highly polar and is typically used as the non-reactive tracer, some partitioning 
into TCE takes place. This occurrence ultimately affects the TCE saturation calculation. 
(The effect of this occurrence will be further discussed in Chapter 4 of this paper.) 
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CHAPTER 4 

ANALYSIS OF RESULTS 
Moment Analysis 

As described in Pope et al. (1994) the “Method of First Moment” theory is typically used 
to determine NAPL saturation in partitioning interwell tracer tests. The retardation factor, 
Rf, and average NAPL saturation, s„, in the subsurface (or “column” during this 
investigation) are both based on the first moments of the partitioning and non-partitioning 
tracers: 

Ks n t 

os) *r- 1+ i rff- 

t m represents the first moment of the non-partitioning tracer [T], 
t p represents the first moment of the partitioning tracer [T], 

K represents the partition coefficient of the partitioning tracer into the NAPL. 

Since figure (18) demonstrates that methanol does not truly represent a “non-reactive” 
tracer, the retardation factor will not reflect the actual retardation which takes place in the 
column or aquifer. The predicted TCE saturation, in turn, will reflect a lower value than 
that which actually exists; this is demonstrated by the analysis described below. 

The first moments of the theoretical non-reactive tracer, methanol, 2-methyl-2-hexanol, 
and 3-methyl-2-hexanol were calculated during this investigation to compare the effects of 
differing moments on the TCE saturation calculation. 

The /th moment, mi, of the effluent response curve is given by the following equation: 

(34) 

o ^ c o y 

For this investigation, the 0 th moment represents the total amount of mass which exits the 
column for the case in which the injected mass equals At pore volumes. The 1 st moment 
represents the time that the center of mass exits the column. 
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For an injection time of 0.5 pore volumes, the theoretical non-reactive tracer should result 
in a 0 th moment of .5 pore volumes and a 1 st moment of 1.25 pore volumes. The 
trapezoidal rule was used to perform the effluent curve integration; thus some error results 
from this approximation. Table VII represents a compilation of each tracer’s 0 th and 1 st 
moments and their associated errors (when they could be quantified). 


Table VII. Tracer/TCE Effluent Moment Approximations 


Tracer 

Moment # 

Value 

% Error 

Non-Reactive 

Mo 

0.5 

0% 


Mj 

1.259999 

0.8% 

Methanol 

Mo 

.500454 

0.0908% 


Mi 

1.811421 

— 

i 2-MethyI-2-Hexanol 

M 0 

.512428 

2.4856% 


M, 

5.418368 

— 

3-Methyl-2-Hexanol 

M 0 

.512498 

2.4996% 


M, 

5.432651 

— 


Effect of Using the “Non-Partitioning” Methanol Assumption 
This portion of the investigation focused on the potential error which results from the 

assumption that methanol acts as a non-partitioning tracer. The method as described 

below was used to quantify the effect on the TCE saturation calculation using methanol as 

the non-reactive tracer. 

First, the retardation factor, Rf, was evaluated using the respective t p and t nr values for 

one partitioning tracer/ideal non-reactive tracer pair. Because all moment calculations 
resulted in systematic moment overestimations, these calculated values were used to 
maintain consistency throughout the error determination. Since the model conditions 
assumed a TCE saturation, s„, of 5% this value was used to calculate a partition 
coefficient, K. The partition coefficient is constant under the linear partitioning 
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assumption and thus can be used in subsequent calculations. Next, the retardation factor 
for each partitioning tracer/methanol pair was calculated. The corresponding TCE 
saturation was then determined by utilizing the partition coefficient found with the 
theoretical tracer and the retardation factor for the tracer/methanol comparison. Finally, 
the error between the “simulated” TCE saturation and that found with methanol was 
calculated. Tables VIII and IX present the results. 


Table VIIL NAPL Saturation Error When Using Methanol as the “Non- 
Reactive” Tracer with 2-Methyl-2-Hexanol 


Ideal 

5.418368 

1.259999 

62.7 

4.30 

.05 

Methanol 

5.418368 

1.811421 

62.7 

2 991 1 

" ' I 

0303119 



•' .• •• ! 


s n ERROR j 

39.38% 


Table IX. NAPL Saturation Error When Using Methanol as the 
“Non-Reactive” Tracer with 3-Methyl-2-Hexanol 



Ideal 

Methanol 




5.432651 I 1.259999 81.9 4.31 

5.432651 | 1.811421 81.9 2.999 

s„ ERROR 


.05 

.0353263 

29.35% 


These results demonstrate that the TCE saturation can be underestimated by almost 40% 
under the given model conditions; thus, small differences in the first moments of an ideal 
non-reactive tracer and actual tracers used for non-reactive tracers may possibly result in 
significant errors. This area deserves further investigation, however, since methanol is 
soluble in water and may not completely follow the investigation assumptions. 
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Effect of Linear Partitioning Assumption 

Another assumption which is typically used in interwell partitioning tracer tests is that the 

tracer/NAPL pair interacts linearly. To analyze this assumption, the Freundlich 

approximation to the UNIFAC-calculated activities was compared to a linear 

approximation to the UNIF AC-calculated activities. The linear approximation was 

x 

determined by solving the equation B’c<f 77^ for B\ This equation is similar to 

equation (15); the only difference is that the Freundlich exponent, A, is set to 1.0. The 
same Co which was determined for each Freundlich run was used for consistency with the 
previously discussed observations. Again, B’ was converted to B using B= B’£. The same 
S values ( see Table v ) were used to convert B’ to B as were used when calculating the 
Freundlich coefficient. The linear parameters determined using the method described 
above are printed in Table X. 


Table X. Linear Approximations to UNIFAC-Calculated 


Tracer/TCE Interactions 


Tracer/TCE Pair 

' 'A 7 

B’ |p 

B 

c« I 

Methanol 

1.0 

.09405 

.55136 

.27262 i 

2-Methyl-2-Hexanol 

1.0 

.24233 

5.14963 

.10581 [ 

3-Methyl-2-Hexanol 

1.0 

.24254 

5.15404 

.10572 j 


Figures (19), (20), and (21) display the results of each tracer when compared with its 
linear approximation. Each comparison also includes a plot of the theoretical “non¬ 
reactive” tracer. 






c/co c/co 
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COMPARISON BETWEEN FREUNDUCH AND LINEAR APPROXIMATIONS FOR 

METHANOL/TCE 

FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 



Figure 19, 


FREUNDUCH VS. LINEAR APPROXIMATIONS FOR 2-METHYL-2-HEXANOL/TCE 
FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 



Figure 20. 
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FREUNDLICH VS. UNEAR APPROXIMATIONS FOR 3-METHYL-2-HEXANOL 
FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 



Figure 21. 


Figures (19), (20), and (21) show that the linear approximations do not yield the same 
effluent response curves as the Freundlich approximations. The 1 st moment calculations 
demonstrate the differences between the time of the center of mass between each tracer 
comparison (see Table XI). 


Table XI: 1 st Moment Differences Between Freundlich 
and Linear Approximations 


Tracer 

Freundlich Mi 

Linear Mi 

AM,. 

Methanol 

1.811421 

1.816846 

.005425 

.■TFP-a:>v+i-..• -••• i 

2-Methyl-2-Hexanol 

5.418368 

6.359433 

.941065 

3-MethyI-2-Hexanol 

5.432651 

6.362715 

**•*♦+♦ v.v-v i v v-v ♦+ H + v « 

.930064 


Figures (19), (20), and (21) simply demonstrate that a linear approximation does not yield 
the same effluent curve as the Freundlich approximation. These results, however, cannot 
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be used to determine the “linear assumption” effects on TCE saturation calculations. To 
correctly quantify the error between TCE saturation calculations, the method as described 
below was used. 

First, each Freundlich effluent curve was treated as the “simulated” column test result. 
Therefore, a linear isotherm must be deduced which will provide the same first moment 
(Mi or t p ) as that produced by the simulated (Freundlich) effluent curve. The 

corresponding linear coefficient, B, was determined by subtracting the first moment of the 
non-reactive tracer from the first moment of the tracer effluent response curve. 
Mathematically, this is represented by the following equation: 

(35) Bljnear (Freundlich) - ^! 1 (non-reactivc)■ 

A corresponding £, value was then obtained using the B obtained by equation (35) and B’ 
obtained by the original linear spreadsheet solver solution. Finally, a corresponding NAPL 
saturation, s„, was determined by solving the following equation: 


(36) s„= 


M.W. T Tjp n +^K / f.W. n ' 


This equation is a combination of equations (6) through (9) and is expressed to solve for 
s„. Table XII presents the difference in TCE saturation determinations: 


Table XII. Linear Assumption Effects on TCE Saturation 
Calculations For Tracer/TCE Pairs 

(a) Methanol: 


Simulated 




(Equivalent Linear 


.82911 i 1.25733 .14143 .0058622 .05 : 

.56142 | 1.0 .094054 .0059694 .04845 

| s n Error: 3.10% ! 
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Table XU (cont.). Linear Assumption Effects on TCE Saturation 
Calculations For Tracer/NAPL Pairs 
(b) 2-Methyl-2-Hexanol: 

* rtit 11 *' * * • '•***-»»* AM ±**M*+ t i* 

il 



[Simulated 
I Equivalent Linear 


9.98927 

4.16837 








.0212505 

.172012 


s„ Error: 22.20% 


(c) 3-Methyl-2-Hexanol 



Simulated 


Equivalent Linear 


*±-. 


9.94644 

4.18265 


1.28392 .46806 .0212505 

1.0 .24254 .0172452 

;:;^:::siSa'Errol*t|P 


03899 

22 , 02 % 


Figures (22), (23), and (24) compare the “simulated” (Freundlich) effluent curves to the 
linear approximations for the tracer/TCE pairs. In addition, comparisons between the 1 st 
moments are also shown. In each comparison, the 1 st moment of the linear 
approximations do not match that of the Freundlich simulated effluent curve exactly. 

These disparities can be attributed to the fact that the calculations utilized the trapezoidal 
approximation to determine the 1 st moments. The errors between the two approximations, 
however, are all within 0.5 %. Thus, the s„ comparisons between the tracer effluent 
curves compared with their linear approximations represent the expected errors which 
result from a linear analysis of partitioning tracer data. 






























LINEAR APPROXIMATION TO FREUNDLICH 1ST MOMENT FOR METHANOUTCE SYSTEM 



TIME (PORE VOLUMES) 


LIN C/CO 
FREUNDLICH 


Figure 22. 


LINEAR APPROXIMATION TO FREUNDLICH 1ST MOMENT FOR 2-METHYL-2-HEXANOL/TCE 

SYSTEM 



LIN C/CO 
FREUNDLICH 


Figure 23. 
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LINEAR APPROXIMATION TO FREUNDLICH 1ST MOMENT FOR 3-METHYL-2-HEXANOL/TCE 

SYSTEM 



Table XII combined with figures (22), (23), and (24) show that a linear approximation to 
non-linear tracer/NAPL interactions can lead to the systematic underestimation of oil 
phase saturations. For the tracers analyzed in this study, the error is as great as 22% when 
the “given” NAPL phase saturation (s„) is 5%; this error applies to tracer mole fractions 
between 0 and 0.05. For larger mole fractions (i.e., x= 0 to 0.20), this error may be larger 
since the system is more non-linear at higher tracer/NAPL mole fractions. These findings 
suggest that a linear approximation to tracer/NAPL interactions may not be appropriate to 
quantify the volume of a NAPL in the subsurface. 

Sensitivity Analysis 

Several variables in the Freundlich approximation to the tracer/TCE interactions were 
modified to determine the sensitivity of changes to these variables on the effluent response 
curve. The tracer 2-methyl-2-hexanol was used as the partitioning tracer in all sensitivity 
analysis tests. First, the Freundlich exponent, A, was changed from its original value of 
1.28623 to a value of 2.0 and then to a value of 3.0. For each exponent, the same 
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injection concentration,^, and the same Freundlich coefficient, B, were used as determined 
by the original solution. Figure (25) displays the effects of these changes. 


FREUNDLICH EXPONENT COMPARISON FOR 2-METHYL-2-HEXANOL/TCE SYSTEM 
FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 



-NON-REACTIVE 

-A=1.28623 

!-A=2.0 

j-A=3.0 __ 


Figure 25. 


The results show that as the exponent, A, is increased, the effluent response curve 
becomes less retarded. A mathematical explanation for this result emanates from the 
characteristic solution for the unfavorable Freundlich isotherm. Equation (21) given in 
the “Advective Transport Equation” section represents the slope of an isochore in the T-t 
plane. The slope of each characteristic curve in the t-t plane is altered when the 

dt 

Freundlich exponent, A, changes. For A>1 and Co concentrations less than 1, — (which 

represents retardation) decreases as A increases. Thus, the behavior of a more non-linear 
system simulates the behavior of a non-reactive tracer; these non-linear systems may not 
be able to correctly predict NAPL saturation since the first moment of the effluent curve 
decreases as the exponent increases. 
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The second sensitivity analysis was conducted by varying the Freundlich coefficient, B. 
For 2-methyl-2-hexanol, the original Freundlich value of B was 9.98927. For the 
sensitivity analysis, B was reduced to one-half its original value. The results are shown in 
figure (26). 


FREUNDLICH COEFF. COMPARISON FOR 2-METHYL-2-HEXANOL/TCE SYSTEM 
FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 



Figure 26. 


Since higher values of the Freundlich coefficient, B, correspond to higher TCE amounts 
present in the column (or subsurface), the effluent results display the expected behavior; 
the 2-methyl-2-hexanol experiences more retardation. In addition, more “smearing” of the 
tracer occurs with higher TCE amounts in the column. These results can also be 
mathematically demonstrated by the characteristic solution to the unfavorable Freundlich 

dt 

isotherm. As the value of B decreases, — also decreases. 

dr 


Finally, the effect of varying Co was analyzed using the original Freundlich parameters for 
the 2-methyl-2-hexanol/TCE system. The Co values which correspond to tracer mole 
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fractions of 0.025 and 0.05 were used in the COMB ALL runs. Figure (27) displays the 
results. 


VARYING MOLE FRACTION COMPARISON, 2-METHYL-2-HEXANOUTCE SYSTEM 
FRACTIONAL BREAKTHROUGH CURVE GENERATED BY THE A-D EQUATION 



Figure 27. 


As can be seen by figure (27), the effluent response curve for the higher mole fraction is 
retarded more than that for the base case. This is expected; if a higher Co is injected into 
the same system for the same time interval, the slope of the system’s characteristic curves 
will be greater. These results are important in field applications since the slightly lower 
peak and the wider base for the higher Co indicates that more “smearing’ occurs when the 
tracer mole fraction is increased. Thus, careful consideration must be taken in the Co 
determination during field applications; although a higher Co allows for analytic instrument 
detection of the effluent, the higher Co may result in the manifestation of more non-linear 


effects. 
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CHAPTER 5 

LANGMUIR APPROXIMATION 


Langmuir isotherms were also fit to the tracer/TCE UNIFAC curves. The “classic’ 
Langmuir isotherm is often expressed in terms of sorbed concentrations: 


(37) 


q represents the sorbed concentration (mass sorbed/solution volume) [ML ]. 
c represents the aqueous phase concentration (mass in solution/solution volume) 
[ML' 3 ]. 


Q represents the adsorption capacity (mass sorbed/solution volume) [ML" 3 ]. 

K represents the selectivity coefficient (mass in solution/solution volume) [ML 3 ]. 


A Langmuir type approximation was initially modeled using the UNIFAC data with 
equation (38): 


(38) x = 


Qc 


K'+c 

x represents the tracer mole fraction 


moles tracer 


moles_tracer + moles_ TCE 
c represents the aqueous phase concentration (mass in solution/solution volume) 
[ML' 3 ]. 

Q’ and K’ represent the Langmuir-type parameters which must be converted to Q 
and K in the “classic” Langmuir expression. 


The Langmuir parameters Q’ and K’ were determined by minimizing the error between the 
UNIFAC calculated “x-c” relationship and the Langmuir approximated “x-c” relationship; 
The Microsoft® Excel Solver was used to vary the Langmuir parameters until the minimal 
error resulted. 
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Langmuir fits were first attempted using mole fraction intervals between 0 and 0.02. As 
can be seen with figures (28), (29), and (30), the Langmuir approximations correlate 
closely with the UNIFAC data. For 2-methyl-2-hexanol and 3-methyl-2-hexanol, the fits 
almost matched exactly. The methanol Langmuir approximation did not fit as closely as 
did the other two tracers; therefore, the mole fraction interval was reduced slightly to x=0 
to 0.13 and is displayed as figure (31). The RMS errors associated with the Langmuir 
approximations are summarized in Table XIII. 


METHANOUTCE UNIFAC Cw COMPARED WITH LANGMUIR APPROXIMATION 



-UNIFAC ACTIVITY (Cw) 

-CALC. LANGMUIR Cw 


Q’= -.07662 
K’= -1.16724 


Figure 28. 
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2-METHYL-2-HEXANOL/TCE UNIFAC Cw COMPARED WITH LANGMUIR APPROXIMATION 



-UNIFAC ACTIVITY (Cw) 

-CALC. LANGMUIR Cw 


Q’=-.11485 
K’= -.58743 


Figure 29. 


3-METHYL-2-HEXANOUTCE UNIFAC Cw COMPARED WITH LANGMUIR APPROXIMATION 



-UNIFAC ACTIVITY (Cw) 

-CALC. LANGMUIR Cw 


Q’= -.1159 
K’= -.59154 


Figure 30. 
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METHANOL/TCE UNIFAC Cw COMPARED WITH LANGMUIR APPROXIMATION 



-UN I FAC ACTIVITY (Cw) 

-CALC. LANGMUIR Cw 


METHANOL MOLE 
FRACTION=0->0.13 

Q’= -.08897 
K’ =-1.26619 


Figure 31. 


Table X1IL Langmuir Values for Tracer/TCE Unifac Approximations 


Tracer K’ Q’ 

RMS Error 

Methanol (0->0.2) 

-1.16724 -.076672 

9.5824 x |g 

Methanol (0->0.13) 

-1.26619 -.08897 

4.2197 x 10‘ 3 

2-Methyl-2-Hexanol 

-.58743 -.11485 

8.8322 x IQ- 4 

3-Methyl-2-Hexanol 

-.59154 -.1159 

8.6942 x 1 O’ 4 


* Mole Fraction Interval = (HK).2 Unless Otherwise Indicated 


Overall, the Langmuir approximation to solving for the mole fraction using the aqueous 
phase concentration is more accurate at larger mole fraction intervals. The Langmuir 
approach can analyze larger mole fraction intervals since the RMS errors between the 
UNIFAC data and the Langmuir approximation are on the order of KT 4 for x= 0 to 0.2. 

In comparison, the Freundlich errors for the same interval were on the order of 10' 3 for the 
same interval as previously shown on Table V. 
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It must be reiterated that the values represented in figures (28) through (31) are expressed 
in terms of x; hence, the presence of negative numbers is not unreasonable. The 
conversion to the “typically” expressed parameters which are normally associated with q 
can be accomplished using a similar analysis to that presented for the Freundlich parameter 
conversion. This conversion will not be shown in this paper since Langmuir model 
simulations were not conducted in this investigation. 
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CHAPTER 6 
CONCLUSIONS 

This investigation yields several conclusions which deal with non-linear partitioning effects 
on NAPL characterization via partitioning tracer tests. UNIFAC computer simulations 
show that non-linear partitioning behavior does, in fact, exist between tracer/NAPL 
interactions at compositions likely to be experienced during tracer tests. To study the 
effects of these non-linearities, the Freundlich fit to the UNIFAC data provides reasonable 
approximations at low concentrations; these Freundlich fits include the following two 
relations: (1) tracer mole fraction (x) vs. aqueous phase concentration (c), and (2) tracer 
storage in the NAPL phase (q) vs. aqueous phase concentration (c). 

Several errors may result when predicting NAPL saturation values using the currently 
accepted methods. First, the assumption that the non-reactive tracer is “truly” non¬ 
reactive may be inaccurate. The NAPL saturation errors using methanol as the non¬ 
reactive tracer in this investigation were as high as 39 percent. Additional investigation in 
this area needs to be conducted, however, since methanol is miscible in water. A second 
type of error results when using the assumption of linear tracer/NAPL partitioning. The 
NAPL saturation errors using 2-methyl-2-hexanol/TCE and 3-methyl-2-hexanol/TCE 
were approximately 22 percent under the model conditions in this study. 

Freundlich parameter variations also lead to changes in the tracer effluent response. As 
the Freundlich exponent. A, increases, less retardation and less spreading occurs. To 
conserve mass, the peak of the effluent curve is higher as A increases. Therefore, higher 
exponent values yield effluent curves which tend to resemble non-reactive tracers. As the 
Freundlich coefficient, B, increases, increased retardation and spreading of the effluent 
curve occurs. Thus, the peak of the effluent curve decreases as B increases. A higher 
coefficient value indicates that more oil is present in the aquifer. Finally, higher injection 
concentrations, Co, yield effluent curves with increased retardation and spreading. Thus, 
care must be exercised in field applications; a balance must be found in which Co is high 
enough so that the tracer can be detected in the effluent, but low enough to minimize the 
non-linear effects associated with high Co values. 
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Models other than the Freundlich may also be used to determine the non-linear 
partitioning effects on NAPL saturation calculations. The Langmuir model is a promising 
candidate for further investigations; for the tracer/TCE pairs investigated in this study, the 
Langmuir model displayed excellent correlations to the UNIFAC data at tracer mole 
fraction intervals between 0 and 0.2. 
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APPENDICES 


Sample UNIFAC Input and Output Files for 2-Methyl-2-Hexanol/TCE 

a and P Values Determined By Microsoft ® Excel Spreadsheet Solver Solution for 
2-Methyl-2-Hexanol/TCE (With RMS Error Calculations) 

A and B Values Determined By Microsoft ® Excel Spreadsheet Solver Solution for 
2-Methyl-2-Hexanol/TCE (With RMS Error Calculations) 

COMBALL Program 

Sample COMBALL Input and Output File for 2-Methyl-2-Hexanol/TCE System 
Using A and B Values and A-D Fractional Breakthrough Curve 



Appendix A 

Sample UNIFAC Input File for 2-Methyl-2-Hexanol/TCE 
For Tracer Mole Fractions From 0-> 1.0 


'2-methyl-2-hexanol/trichloroethylene ' / <"Text Line 

110 / <-Use “VLE” internal file parameters and write to both screen &“gamma.out” 

2 / <-# of components being compared 

0 / 4-Retrieve “R” and “Q” values from internal data files 

3132141 15/ 4-1“ compound’s# in subgroup, subgroup#... (There are 4 subgroups) 

1 8 3 70 / <-2 nd compound’s# in subgroup, subgroup#... (There are 2 subgroups) 

0 / 

298. / 4-Temperature (K) 

0.0 1.0/ 4-Mole fraction of 1“ compound followed by mole fraction of 2 compound 

298./ 

0.1 .90/ 

.2 .80 / 

298. / 

.3 .70/ 

298 / 

.4 .60 / 

298. / 

.5 .50 / 

298./ 

.6 .40 / 

298. / 

.7 .30 / 

298. / 

.8 .20 / 

298. / 

.9 .10/ 

298. / 

1.0 0 . 0 / 

0.0 


/ ^-Signifies end of data file 


1 


IOUT = 1 

F IOUT = 0 : OUTPUT ONLY ON FILE 
F IOUT = 1 : OUTPUT ON BOTH FILE AND SCREEN 


2-methyl-2-hexanol/trichloroethylene 


MODEL USED FOR LIQUID PHASE: ORIGINAL UNFAC (VLE) 


GROUP SPECIFICATION: 


SUBGROUP MAIN R Q PRESENCE IN COMPONENT NO. 
NO. NAME GROUP 1 2 

1 CH3. 1 0.9011 0.8480 3 0 

2 CH2. 1 0.6744 0.5400 3 0 

4 C. 1 0.2195 0.0000 1 0 

8 CHC. 2 0.8886 0.6760 0 1 

15 OH. 5 1.0000 1.2000 1 0 

70 CL-(C=C) 37 0.7910 0.7240 0 3 


PARAMETER COEFFICIENTS FOR INTERACTIONS 
1 2 5 37 

1 


0.0000 

86.02 

986.5 

-4.189 

-35.36 

0.0000 

524.1 

-66.46 

156.4 

457.0 

0.0000 

225.8 

47.41 

124.2 

738.9 

0.0000 


MOLECULAR PARAMETERS 
COMP. R Q 

1 5.9460 5.3640 

2 3.2616 2.8480 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
0.000000 5.14722 








100.000000 


1.00000 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
10.000000 2.72672 

90.000000 1.03114 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
20.000000 1.87894 

80.000000 1.09943 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
30.000000 1.48710 

70.000000 1.18730 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
40.000000 1.27762 

60.000000 1.28741 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
50.000000 1.15647 

50.000000 1.39583 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
60.000000 1.08387 

40.000000 1.51006 


TEMPERATURE: 298.000 



MOLE PER CENT ACTIVITY COEFFICIENT 
70.000000 1.04043 

30.000000 1.62835 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
80.000000 1.01567 

20.000000 1.74932 


TEMPERATURE: 298.000 

MOLE PER CENT ACTIVITY COEFFICIENT 
90.000000 1.00346 

10.000000 1.87190 


TEMPERATURE: 298.000 


MOLE PER CENT ACTIVITY COEFFICIENT 
100.000000 1.00000 
0.000000 1.99518 



Appendix B 

Excel Spreadsheet Solver Solution For 2-Methyl-2-Hexanol/TCE 
(With RMS Error Calculation) 


2-METHYL-2-HEXANOL/TCE UNIFAC GENERATED ACTIVITIES COMPARED WITH 
FREUNDLICH APPROXIMATIONS, "BETA" AND "ALPHA" OBTAINED BY 

EXCEL SOLVER ROUTINE ___ 

BETA= 0.41835494 

ALPHA= 1.24654054 


FREUNDLICH APPROXIMATION 


UNIFAC 

MOLE 

CALC. 

UNIFAC FREUND. 


ACTIVITY 

FRACTION 

ACTIVITY MOLE 


(Cw) 

TRACER 

(Cw) FRACTION 

ERROR 

0 

0 

0 0 

0 

0.005103932 

0.001 

0.005103932 0.00058124 

1.7536E-07 

0.01012258 

0.002 

0.01012258 0.00136477 

4.03517E-07 

0.015057834 

0.003 

0.015057834 0.002238989 

5.79138E-07 

0.019911536 

0.004 

0.019911536 0.003171831 

6.85864E-07 

0.024685475 

0.005 

0.024685475 0.004146274 

7.28847E-07 

0.029381406 

0.006 

0.029381406 0.005151518 

7.19921E-07 

0.034001023 

0.007 

0.034001023 0.006180023 

6.72361 E-07 

0.038545984 

0.008 

0.038545984 0.00722621 

5.98751 E-07 

0.043017912 

0.009 

0.043017912 0.00828578 

5.101 IE-07 

0.04741838 

0.01 

0.04741838 0.009355323 

4.15608E-07 

0.051748917 

0.011 

0.051748917 0.010432075 

3.22539E-07 

0.05601102 

0.012 

0.05601102 0.011513757 

2.36432E-07 

0.060206146 

0.013 

0.060206146 0.012598467 

1.61228E-07 

0.064335726 

0.014 

0.064335726 0.013684605 

9.94739E-08 

0.068401125 

0.015 

0.068401125 0.0147708 

5.25328E-08 

0.072403728 

0.016 

0.072403728 0.015855892 

2.07672E-08 

0.076344824 

0.017 

0.076344824 0.016938868 

3.73715E-09 

0.08022573 

0.018 

0.08022573 0.018018867 

3.55966E-10 

0.084047659 

0.019 

0.084047659 0.019095124 

9.04862E-09 

0.08781186 

0.02 

0.08781186 0.020166991 

2.7886E-08 

0.091519533 

0.021 

0.091519533 0.021233899 

5.47087E-08 

0.095171824 

0.022 

0.095171824 0.022295345 

8.72288E-08 

0.098769889 

0.023 

0.098769889 0.023350902 

1.23132E-07 

0.102314832 

0.024 

0.102314832 0.024400191 

1.60153E-07 

0.105807725 

0.025 

0.105807725 0.02544288 

1.96142E-07 

0.109249634 

0.026 

0.109249634 0.026478686 

2.29141 E-07 

0.112641597 

0.027 

0.112641597 0.027507366 

2.57421 E-07 

0.115984624 

0.028 

0.115984624 0.028528707 

2.79531 E-07 

0.119279668 

0.029 

0.119279668 0.029542518 

2.94326E-07 

0.12252774 

0.03 

0.12252774 0.030548659 

3.01026E-07 

0.125729707 

0.031 

0.125729707 0.031546976 

2.99183E-07 

0.12888656 

0.032 

0.12888656 0.032537385 

2.88782E-07 

0.131999142 

0.033 

0.131999142 0.033519778 

2.70169E-07 

0.135068366 

0.034 

0.135068366 0.034494096 

2.44131 E-07 

0.138095055 

0.035 

0.138095055 0.035460273 

2.11852E-07 

0.141080076 

0.036 

0.141080076 0.036418277 

1.74956E-07 

0.144024202 

0.037 

0.144024202 0.037368065 

1.35472E-07 




0.146928292 

0.149793033 

0.15261928 

0.155407712 

0.158159064 

0.160874094 

0.163553456 

0.166197825 

0.168807902 

0.171384325 

0.173927664 

0.176438661 

0.1789178 


Appendix B 

Excel Spreadsheet Solver Solution For 2-Methyl-2-Hexanol/TCE 
(With RMS Error Calculation) 


0.038 

0.146928292 

0.038309639 

9.5876E-08 

0.039 

0.149793033 

0.039242961 

5.90303E-08 

0.04 

0.15261928 

0.040168066 

2.82462E-08 

0.041 

0.155407712 

0.041084942 

7.21514E-09 

0.042 

0.158159064 

0.04199361 

4.08313E-11 

0.043 

0.160874094 

0.042894111 

1.12124E-08 

0.044 

0.163553456 

0.043786464 

4.55978E-08 

0.045 

0.166197825 

0.044670703 

1.08437E-07 

0.046 

0.168807902 

0.045546884 

2.05314E-07 

0.047 

0.171384325 

0.04641505 

3.42167E-07 

0.048 

0.173927664 

0.04727523 

5.25292E-07 

0.049 

0.176438661 

0.048127519 

7.61224E-07 

0.05 

0.1789178 

0.048971934 

1.05692E-06 


IERROR: 1.32774E-051 

RMS= 0.000515314 


EQUATION USED FOR RMS ERROR-> 


(x-xf 


RMS ERROR = 


'# DATA POINTS 






Appendix C 

Excel Spreadsheet Solver Solution for 2-Methyl-2-Hexanol 
(With RMS Error Calculation) 


2-METHYL-2-HEXANOL/TCE UNIFAC GENERATED ACTIVITIES COMPARED WITH 
FREUNDLICH APPROXIMATIONS, "B" AND "A" OBTAINED BY EXCEL SOLVER ROUTINE 


BETA= 

ALPHA= 

0.418355 

1.246541 


fit 0.47007241 

A- 1.28623477 

ZETA (g/i)= 
B 

21.2505 

9.989273676 

UNIFAC 

ACTIVITY 

(Cw) 

MOLE 

FRACTION 

TRACER 

X/(1-X) 

FREUNDLICH APPROXIMATION 

CALC. X/(1- 
X) 

ERROR 

0 

0 

0 

0 


0 

0.005103932 

0.001 

0.001001001 

0.000529654 


2.22168E-07 

0.01012258 

0.002 

0.002004008 

0.001277912 


5.27215E-07 

0.015057834 

0.003 

0.003009027 

0.002129803 


7.73034E-07 

0.019911536 

0.004 

0.004016064 

0.003050803 


9.31729E-07 

0.024685475 

0.005 

0.005025126 

0.004022232 


1.0058E-06 

0.029381406 

0.006 

0.006036217 

0.005032067 


1.00832E-06 

0.034001023 

0.007 

0.007049345 

0.006071817 


9.55561 E-07 

0.038545984 

0.008 

0.008064516 

0.007135131 


8.63757E-07 

0.043017912 

0.009 

0.009081736 

0.008217071 


7.47646E-07 

0.04741838 

0.01 

0.01010101 

0.009313681 


6.19886E-07 

0.051748917 

0.011 

0.011122346 

0.010421731 


4.90861 E-07 

0.05601102 

0.012 

0.012145749 

0.011538533 


3.68712E-07 

0.060206146 

0.013 

0.013171226 

0.012661826 


2.59488E-07 

0.064335726 

0.014 

0.014198783 

0.013789692 


1.67356E-07 

0.068401125 

0.015 

0.015228426 

0.014920473 


9.48351 E-08 

0.072403728 

0.016 

0.016260163 

0.016052756 


4.30173E-08 

0.076344824 

0.017 

0.017293998 

0.017185297 


1.1816E-08 

0.08022573 

0.018 

0.018329939 

0.018317024 


1.66796E-10 

0.084047659 

0.019 

0.019367992 

0.019446982 


6.23948E-09 

0.08781186 

0.02 

0.020408163 

0.02057435 


2.7618E-08 

0.091519533 

0.021 

0.02145046 

0.021698399 


6.1474E-08 

0.095171824 

0.022 

0.022494888 

0.022818481 


1.04713E-07 

0.098769889 

0.023 

0.023541453 

0.023934035 


1.5412E-07 

0.102314832 

0.024 

0.024590164 

0.025044557 


2.06473E-07 

0.105807725 

0.025 

0.025641026 

0.026149602 


2.5865E-07 

0.109249634 

0.026 

0.026694045 

0.027248783 


3.07734E-07 

0.112641597 

0.027 

0.027749229 

0.028341758 


3.5109E-07 

0.115984624 

0.028 

0.028806584 

0.029428224 


3.86435E-07 

0.119279668 

0.029 

0.029866117 

0.030507905 


4.11891E-07 

0.12252774 

0.03 

0.030927835 

0.031580585 


4.26083E-07 

0.125729707 

0.031 

0.031991744 

0.032646038 


4.281 E-07 

0.12888656 

0.032 

0.033057851 

0.033704111 


4.17652E-07 

0.131999142 

0.033 

0.034126163 

0.034754636 


3.94977E-07 

0.135068366 

0.034 

0.035196687 

0.035797494 


3.60969E-07 

0.138095055 

0.035 

0.03626943 

0.036832565 


3.17121 E-07 

0.141080076 

0.036 

0.037344398 

0.037859768 


2.65606E-07 

0.144024202 

0.037 

0.038421599 

0.038879011 


2.09226E-07 

0.146928292 

0.038 

0.03950104 

0.039890254 


1.51488E-07 




0.149793033 

0.15261928 

0.155407712 

0.158159064 

0.160874094 

0.163553456 

0.166197825 

0.168807902 

0.171384325 

0.173927664 

0.176438661 

0.1789178 


Appendix C 

Excel Spreadsheet Solver Solution for 2-Methyl-2-Hexanol 


0.039 

(With RMS Error Calculation) 
0.040582726 0.040893417 

9.65289E-08 

0.04 

0.041666667 

0.041888498 

4.92089E-08 

0.041 

0.042752868 

0.042875447 

1.50257E-08 

0.042 

0.043841336 

0.043854254 

1.66862E-10 

0.043 

0.044932079 

0.044824929 

1.14813E-08 

0.044 

0.046025105 

0.04578746 

5.64748E-08 

0.045 

0.047120419 

0.046741857 

1.43309E-07 

0.046 

0.048218029 

0.04768815 

2.80772E-07 

0.047 

0.049317943 

0.048626359 

4.78289E-07 

0.048 

0.050420168 

0.049556488 

7.45943E-07 

0.049 

0.051524711 

0.050478617 

1.09431E-06 

0.05 

0.052631579 

0.051392739 

1.53472E-06 


IERROR: 1.88453E-051 

RMS= 0.000613926 


EQUATION USED FOR RMS ERROR: 


RMS ERROR = 







Appendix D 
COMBALL Program 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


********+********++*****+***************************************** 
* * 

* THIS PROGRAM DETERMINES THE EFFLUENT CURVE FOR A MULTI- * 

* WELL INJECTION TEST. TO RUN THE MODEL IT IS NECESSARY TO * 

* ASSUME AN ISOTHERM, AND ITS PARAMETERS, AS WELL AS THE FRAC- * 

* TIONAL BREAKTHROUGH CURVE. * 

* COMBALL WAS WRITTEN BY WILLIAM R. WISE AT THE UNIVERSITY * 

* OF TEXAS AT AUSTIN IN THE SUMMER OF 1984, BASED ON THE THEORY * 

* OF THE COMBINATION DIAGRAM ALSO DEVELOPED BY THE AUTHOR. * 

* * 
****************************************************************** 




* 

★ 

★ 

DATA ENTRY INSTRUCTIONS 


★ 

+ 

* 

* VARIABLE 

DEFINITION 

FORMAT 

★ 

★ 

* 

* 

* ISOTHER 

ISOTHERM (EITHER LANGMUIR OR 

A10 

★ 

★ 

+ 

FREUNDLICH) 


-*■ 

★ 

* K OR B 

K IF ISOTHER=LANGMUIR 

F15.5 

* 

★ 

B IF ISOTHER=FREUNDLICH 


* 

* 

* Q OR A 

Q IF ISOTHER=LANGMUIR 

F15.5 

* 

* 

A IF ISOTHER=FREUNDLICH 


+ 

* 

* CO 

INJECTED CONCENTRATION 

F15.5 

* 

★ 

* DELT 

INJECTION TIME 

F15.5 

■k 

* 

* TMAX 

MAXIMUM TIME OF EVALUATION 

F15.5 

* 

* 

* TINCR 
£ 

TIME INCREMENT OF EVALUATION 

F15.5 

* 

* 

* NDIV 

NUMBER OF DIVISIONS FOR INTEGRATING 

15 

* 

+ 

ACROSS SPREADING WAVE 


★ 

+ 

* NF 

NUMBER OF COORDINATES OF FRACTIONAL 

15 

* 

★ 

BREAKTHROUGH CURVE 


* 

★ 

* TFPOINT(J), 

JTH COORDINATE OF F CURVE 

2F15.5 

* 

* FPOINT (J) 



★ 


****************************************************************** 


PROGRAM COMBALL 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 
C TFPOINT(200),FPOINT(200), 

C CE(200) , PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL K 

CHARACTER*10 ISOTHER,INPUT,OUTDAT 

WRITE (*,*) 'PLEASE ENTER INPUT DATA FILE NAME* 

READ (*,*) INPUT 

WRITE (*,*) 'PLEASE ENTER OUTPUT DATA FILE NAME' 

READ (*,*) OUTDAT 
OPEN (1,FILE=INPUT) 

OPEN (2,FILE=OUTDAT) 

CALL ASSEMBL 

ITMAX = IFIX(TMAX / TINCR) 

IF (ISOTHER.EQ.'FREUNDLICH') THEN 
IF (A.LT.l) THEN 

DO 110 J=l,ITMAX 
T = J * TINCR 







CE(J) = CEFOl(T) 

110 CONTINUE 

ELSE IF (A.EQ.l) THEN 
DO 120 J=l,ITMAX 
T = J * TINCR 
CE(J) = CEF1(T) 

120 CONTINUE 

ELSE 

DO 130 J=l,ITMAX 
T = J * TINCR 
CE(J) = CEFIUP(T) 

WRITE(2,125) T,CE(J) 

125 FORMAT(2F15.5) 

130 CONTINUE 

END IF 
ELSE 

DO 140 J=l,ITMAX 
T = J * TINCR 
CE (J) = CEL (T) 

140 CONTINUE 

END IF 

CALL MASSBAL(CE,TINCR,DELT,ITMAX,CMASSIN,CMASSOU,PERCENT,MASSREM, 

C CO) 

CALL OUTPUT 
END 
C 

c ****************************************************************** 

C * * 

c * THE ASSEMBL SUBROUTINE READS THE INPUT DATA. * 

C * * 

C ****************************************************************** 

c 

SUBROUTINE ASSEMBL 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 

C TFPOINT(200),FPOINT(200), 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL K 

CHARACTER*10 ISOTHER 
READ (1,500) ISOTHER 
500 FORMAT (A10) 

IF (ISOTHER.EQ.’FREUNDLICH’) THEN 

READ (1,510) B, A, CO,DELT 

ELSE IF (ISOTHER.EQ.'LANGMUIR ') THEN 

READ (1,510) K,Q,C0,DELT 

ELSE 

WRITE (2,515) ISOTHER 

515 FORMAT ('1',///,20X,'EXECUTION WAS TERMINATED',///, 

C 20X,A10,' IS NOT AN ACCEPTABLE ISOTHERM',///, 

C 2OX,'ONLY LANGMUIR OR FREUNDLICH MAY BE USED') 

STOP 
END IF 

READ (1,510) TMAX 
READ (1,510) TINCR 
510 FORMAT (F15.5) 

READ (1,520) NDIV 
READ (1,520) NF 
520 FORMAT (15) 

DO 540 J=l,NF 

READ (1,530) TFPOINT(J), FPOINT(J) 

530 FORMAT (2F15.5) 

WRITE(2,535) TFPOINT(J),FPOINT(J) 

535 FORMAT(2F15.5) 

540 CONTINUE 

RETURN 
END 
C 

C ****************************************************************** 

c * * 

C * THE CEL FUNCTION SUBROUTINE DETERMINES THE EFFLUENT * 

C * CONCENTRATION AT THE PRODUCTION WELL AT A TIME T FOR THE * 

C * LANGMUIR ISOTHERM. * 

C * * 




c 

c 




FUNCTION CEL(T) 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 

C TFPOINT(200),FPOINT(200), • 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL K 

CHARACTER*10 ISOTHER 
TO = TF(0.) 

T1 = TF(1.) 

IF (T.GT.(1 + Q / (K + CO)) * TO) GO TO 200 
CEL = 0. 

RETURN 

200 IF (T.GT.DELT + (1 + K * Q / (K+ CO)**2) * TO) GO TO 210 
TFSTEP = T / (1 + Q / (K + CO)) 

CEL = CO * F(TFSTEP) 

RETURN 

210 TPLAT = DELT * (K + CO + Q) * (K + CO) / (Q * CO) 

TFO = (T - DELT) / (1 + Q / K) 

IF (TF0.GT.T1) RETURN 
IF (TFO.LT.TO) TFO = TO 
IF (T. LE.TPLAT) THEN 

TFSTOP = (T - DELT) / (1 + K * Q / (K+ C0)**2) 

IF (TFSTOP.GT.T1) TFSTOP = T1 
ELSE 

ALPHA « (1 + Q / K)**2 
BETA = 4 * DELT * CO / K 

C -2* (1+Q/K) * (T- DELT * (1 - CO / K)) 

GAMMA = (T - DELT * (1 - CO / K))**2 
C - 4 * DELT * CO * (T - DELT) / K 

TFSTOP = (-BETA + SQRT(BETA**2 - 4 * ALPHA * GAMMA)) / (2 * ALPHA) 
IF (TFSTOP.GT.T1) TFSTOP = T1 
END IF 

IF (TFSTOP.GT.TO) GO TO 220 
CEL = 0. 

RETURN 

220 TF1 = TFO 
FI = F(TF1) 

DELF = F(TFSTOP) - FI 

C = (SQRT(K * Q * TF1 / (T - DELT - TF1)) - K) / 2 
C = C + (SQRT(K * Q * TFSTOP / (T - DELT - TFSTOP)) - K) / 2 
DO 230 J=2,NDIV 
F2 = FI + DELF / NDIV 
TF2 = TF(F2) 

C = C + SQRT(K * Q * TF2 / (T - DELT - TF2)) - K 
FI = F2 

230 CONTINUE 

CEL = C * DELF / NDIV 
IF (T.GE.TPLAT) RETURN 
TFSTEP = T / (1 + Q / (K + CO)) 

IF (TFSTEP.GT.T1) TFSTEP = T1 

CEL = CEL + CO * (F(TFSTEP) - F (TFSTOP)) 

RETURN 

END 

C 

C ****************************************************************** 

C * * 

C * THE CEF01 FUNCTION DETERMINE THE EFFLUENT CONCENTRATION * 

C * FOR THE FREUNDLICH ISOTHERM WHEN THE EXPONENT IS LESS THAN 1. * 

C * * 

c ****************************************************************** 

c 

FUNCTION CEFOl(T) 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 

C TFPOINT(200),FPOINT(200), 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL IOTA 

CHARACTER*10 ISOTHER 
TO = TF(0.) 

T1 = TF(1.) 

ZETA = A - 1 
ETA = B * CO**ZETA 



nnnoooon 


IF (T.-4E. (1 + ETA) * TO) THEN 
CEF01 = 0. 

ELSE IF (T.LE.DELT + (1 + A * ETA) * TO) THEN 
TFSTEP = T / (1 + ETA) 

IF (TFSTEP.GT.T1) TFSTEP = T1 
CEF01 = CO * F(TFSTEP) 

ELSE 

TPLAT = (1 + ETA) * DELT / (ETA * (-ZETA)) 

IF (T.LE.TPLAT) THEN 

TFSTOP = (T-DELT) / (1 + A * ETA) 

IF (TFSTOP.GT.T1) TFSTOP = T1 
ELSE 

R = T / (1 + ETA) 

IOTA = DELT * CO * (A**(A/ZETA) ) * (B**(1/ZETA)) / (-ZETA) 

DO 310 J=l,1000 

G = R**(ZETA/A) - (T - DELT) / R**(l/A) + IOTA**(ZETA/A) 

IF (ABS(G).LE.10E-8) GO TO 320 

DGDR = (ZETA / A) / R** (1/A) + (1 / A) * (T - DELT) / 

C R**((A+1)/A) 

R = R - G / DGDR 
310 CONTINUE 

STOP 

320 TFSTOP = R 

IF (TFSTOP.GT.T1) TFSTOP = T1 
END IF 

IF (TFSTOP.LE.TO) THEN 
CEF01 = 0. 

ELSE 
TF1 = TO 
FI = F(TF1) 

DELF = F(TFSTOP) - FI 

C = (((T - DELT - TF1) / (A * B * TF1))**(1/ZETA)) / 2 
C = C + (((T - DELT - TFSTOP) / (A * B * TFSTOP))**(1/ZETA) ) / 2 
DO 330 J = 2,NDIV 
F2 = FI + DELF / NDIV 
TF2 = TF(F2) 

C = C + ((T - DELT - TF2) / (A * B * TF2))**(1/ZETA) 

FI = F2 

330 CONTINUE 

CEF01 = C * DELF / NDIV 
IF (T.GE.TPLAT) RETURN 
TFSTEP = T / (1 + ETA) 

IF (TFSTEP.GT.T1) TFSTEP = T1 

CEF01 = CEF01 + CO * (F(TFSTEP) - F(TFSTOP)) 

END IF 
END IF 
RETURN 
END 

****************************************************************** 


* * 

* THE CEF1 FUNCTION DETERMINES THE EFFLUENT CONCENTRATION * 

* FOR THE FREUNDLICH ISOTHERM WHEN THE EXPONENT EQUALS 1. * 

* * 




FUNCTION CEF1(T) 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 
C TFPOINT(200),FPOINT(200) , 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

CHARACTER*10 ISOTHER 
TO = TF(0.) 

T1 = TF(1.) 

TFSTEP = T / (1 + B) 

IF (TFSTEP.GT.T1) TFSTEP = T1 
IF (T.LE.(1 + B) * TO) THEN 
CEF1 = 0. 

ELSE IF (T.LE.DELT + (1 + B) * TO) THEN 
CEF1 = CO * F(TFSTEP) 

ELSE 

TFCO = (T - DELT) / (1 + B) 

IF (TFC0.GT.T1) TFCO = T1 



CEF1 = CO * (F(TFSTEP) - F(TFCO)) 

END IF 
RETURN 
END 
C 

c ****************************************************************** 

c * + 

c * THE CEF1UP FUNCTION DETERMINES THE EFFLUENT CONCENTRATION * 

C * FOR THE FREUNDLICH ISOTHERM WHEN THE EXPONENT EXCEEDS 1. * 

C * * 

c ***★**★★★*★★*★*★**★★★★*★★★*★★★★**★+★****+★****★★★**★*★★★★**★★*★★** 

c 

FUNCTION CEFIUP(T) 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 

C TFPOINT(200),FPOINT(200), 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL IOTA 

CHARACTER*10 ISOTHER 
TO = TF(0.) 

T1 = TF(1.) 

ZETA = A - 1 
ETA = B * CO**ZETA 
IF (T.LE.TO) THEN 
CEF1UP = 0. 

ELSE 

TFPLAT = DELT / {ETA * ZETA) 

TPLAT = (1 + A * ETA) * TFPLAT 

IF (T.LT.TPLAT) THEN 

IF (TO.LT.TFPLAT) THEN 

IF (T.LT.(1 + A * ETA) * TO) THEN 

TF1 = TO 

ELSE 

TF1 = T / (1 + A * ETA) 

END IF 
ELSE 

TF1 = TO 
END IF 
ELSE 

R = (T - DELT) / (1 + ETA) 

IOTA = DELT * CO * (A** (A/ZETA)) * (B**(1/ZETA)) / ZETA 
DO 410 J=l,100000 

G= (R** (ZETA/A) ) - (T / (R** (1/A) ) ) + (IOTA** (ZETA/A) ) 

IF (ABS(G).LE.10E-6) GO TO 420 

DGDR = ((ZETA/A) / (R**(1/A) )) + (1/A) * T / (R**((A+l)/A)) 

R = R - G / DGDR 
410 CONTINUE 

STOP 

420 TF1 = R 

IF (TF1.GT.T1) THEN 
CEF1UP = 0. 

RETURN 

ELSE 

CONTINUE 
END IF 

IF (TF1.LT.T0) TF1 = TO 
END IF 
TFSTOP = T 

IF (TFSTOP.GT.T1) TFSTOP = T1 
FI = F(TF1) 

DELF = F(TFSTOP) - FI 

C = (((T - TF1) / (A * B * TF1))**(1/ZETA)) / 2 
C = C + (((T - TFSTOP) / (A * B * TFSTOP))**(1/ZETA)) / 2 
DO 430 J=2,NDIV 
F2 = FI + DELF / NDIV 
TF2 = TF(F2) 

C = C + ((T - TF2) / (A * B * TF2))**(1/ZETA) 

FI = F2 

430 CONTINUE 

CEF1UP = C * DELF / NDIV 
IF (T.GE.TPLAT) RETURN 
IF (TO.GE.TFPLAT) RETURN 
IF (T.LE.(1 + A * ETA) * TO) RETURN 



c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


TFSTEP = (T- DELT) / (1 + ETA) 

IF (TFSTEP.LT.TO) TFSTEP = TO 
IF (TFSTEP.GT.T1) TFSTEP = T1 
TFCO = T / (1 + A * ETA) 

IF (TFC0.GT.T1) TFCO = T1 

CEF1UP = CEF1UP + CO * (F(TFCO)- F(TFSTEP)) 

END IF 
RETURN 
END 

**★★★★★**★★*★*★★*★★★*★★*★**★*★*★*★★*******+*********************** 
* * 

* THE F FUNCTION SUBROUTINE DETERMINES THE F VALUE FOR A * 

* GIVEN TF (THE ARGUMENT TH) THROUGH LINEAR INTERPOLATION. * 

* * 
*******************★+*************++********+********************* 


FUNCTION F(TH) 

COMMON ISOTHER,B, A,K, Q,CO,DELT, TMAX,TINCR,NDIV,NF, ITMAX, 

C TFPOINT(200),FPOINT(200), 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL K 

CHARACTER*10 ISOTHER 
IF (TH.GT.TF(0.)) GO TO 700 
F = 0. 

RETURN 
700 J = 2 

710 IF (TH - TFPOINT(J)) 740,730,720 
720 J = J + 1 

IF (J.LT.NF) GO TO 710 
WRITE (2,777) 

777 FORMAT (’1',///,20X,’EXECUTION WAS TERMINATED' , 111 , 2OX, 

C ’FRACTIONAL BREAKTHROUGH CURVE IS INSUFFICIENT',///, 

C 20X,'FOR GENERATED TF ARGUMENT FOR FUNCTION F') 

STOP 

730 F = FPOINT(J) 

RETURN 

740 F = FPOINT(J-l) + ((FPOINT(J) - FPOINT(J-l)) / (TFPOINT(J) - 
C TFPOINT(J-l))) * (TH - TFPOINT(J-l)) 

RETURN 

END 

****************************************************************** 


* THE TF FUNCTION SUBROUTINE DETERMINES THE VALUE OF TF FOR * 

* A CORRESPONDING VALUE OF F (THE ARGUMENT FH) THROUGH LINEAR * 

* INTERPOLATION. * 

* * 

****************************************************************** 


FUNCTION TF(FH) 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 

C TFPOINT(200),FPOINT(200), 

C CE(200) , PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL K 

CHARACTER*10 ISOTHER 
IF (FH.NE.0) GO TO 800 
TF = TFPOINT(1) 

RETURN 
800 J = 2 

IF (FH.EQ.l..AND.FPOINT(NF).LT.l.) THEN 
TF = 10 * TFPOINT(NF) 

RETURN 
ELSE 
END IF 

810 IF (FH - FPOINT(J)) 840,830,820 
820 J = J + 1 

IF (J.LT.NF) GO TO 810 
WRITE (2,888) 

888 FORMAT ('1',///,20X,'EXECUTION WAS TERMINATED',///,2OX, 

C 'FRACTIONAL BREAKTHROUGH CURVE IS INSUFFICIENT’,///, 

C 20X,'FOR GENERATED F ARGUMENT FOR FUNCTION TF') 



c 

c 

c 

c 

c 

c 

c 

c 


STOP 

830 TF « TFPOINT(J) 

RETURN 

840 TF = TFPOINT(J-l) + ((TFPOINT(J) - TFPOINT(J-l)) / (FPOINT(J) 

C - FPOINT(J-l))) * (FH - FPOINT(J-l)) 

RETURN 

END 

****************************************************************** 


* THE MASSEAL SUBROUTINE PERFORMS A MASS BALANCE ON THE * 

* INJECTED AND PRODUCED MASS. * 

* * 

★+*★*********************+*****★**********★*****+****★*★*★****★*** 


SUBROUTINE MAS SEAL (CE, TINCR, DELT, ITMAX, CMASSIN, CMASSOU, 

C PERCENT,MASSREM, CO) 

DIMENSION CE(200),PERCENT(200) 

CMASSIN = CO * DELT 
CMASSOU = 0. 

M = 1 

DO 900 L=l, ITMAX 

IF (L.GT.1.AND.M.EQ.1) THEN 

IF (CE(L-l).EQ.0..AND.CE(L).GT.0) THEN 

PERCENT(1) = L * TINCR 

M “ 2 

ELSE 

END IF 

ELSE 

END IF 

CMASSOU = CMASSOU + CE(L) 

IF (CMASSOU.GE.(CMASSIN * (M - 1))/10..AND.M.GT.1) THEN 

PERCENT(M) * (L - 1) * TINCR 

M = M + 1 

ELSE 

END IF 

900 CONTINUE 

MASSREM = M - 1 
RETURN 
END 


C 

c ****************************************************************** 

c * * 

C * THE OUTPUT SUBROUTINE PRINTS THE INPUT AND OUTPUT RESULTS.* 

C * ' * 

C ****************************************************************** 

c 


SUBROUTINE OUTPUT 

COMMON ISOTHER,B,A,K,Q,CO,DELT,TMAX,TINCR,NDIV,NF,ITMAX, 

C TFPOINT(200),FPOINT(200), 

C CE(200),PERCENT(200),CMASSIN,CMASSOU,MASSREM 

REAL K 

CHARACTER*10 ISOTHER 

IF (ISOTHER.EQ.'LANGMUIR ’) THEN 

WRITE (2,610) 

610 FORMAT (’1’,////,25X,’LANGMUIR MULTIWELL INJECTION SIMULATION’) 
WRITE (2,620) K,Q, CO,DELT 

620 FORMAT (///,40X,'INPUT DATA',///,35X,’LANGMUIR PARAMETERS',//, 

C 24X,'EQUILIBRIUM CONSTANT (K) = ',F12.3,//, 

C 24X,'ADSORPTION CAPACITY (Q) = ',F12.3,//, 

C 24X,'INJECTED CONCENTRATION (CO) = ',F12.3,//, 

C 24X,'INJECTION TIME (DELT) = ',F12.3,////, 

C 41X,'ISOTHERM',//,32X,'C «= (Q) (C) / (K + C) ', 

C /,33X'SORBED') 

ELSE 

WRITE (2,615) 

615 FORMAT ('1’,////,24X,'FREUNDLICH MULTIWELL INJECTION SIMULATION') 
WRITE (2,625) B,A,CO,DELT 

625 FORMAT (///,40X,'INPUT DATA’,///,34X,'FREUNDLICH PARAMETERS',//, 
C 24X,'FREUNDLICH COEFFICIENT (B) = ',F12.3,//, 

C 24X,'FREUNDLICH EXPONENT (A) = ',F12.3,//, 

C 24X,'INJECTED CONCENTRATION (CO) = ',F12.3,//, 



c 24X,'INJECTION TIME (DELT) = ',F12.3,////, 

C 41X, 'ISOTHERM',//,52X,'A' ,/,36X, 'C = (B) (C) ' , 

C /,37X,'SORBED') 

END IF 

WRITE (2,630) 

630 FORMAT (///,35X,'MASS BALANCE RESULTS',//,31X,’PERCENT',15X, 

C ’TIME',/,31X,'REMOVAL',/) 

DO 632 I = 1,MASSREM 
M = 10 * (I - 1) 

WRITE (2,631) M, PERCENT(I) 

631 FORMAT (33X,13,9X,F13.3) 

632 CONTINUE 

WRITE (2,633) CMASSIN,CMASSOU 

633 FORMAT (//,31X,'MASS INJECTED = 'F11.3,//, 

C 3IX,'MASS REMOVED = 'F11.3) 

WRITE (2,635) 

635 FORMAT ('1’,///,30X,'FRACTIONAL BREAKTHROUGH CURVE',///, 

C 34X,’F’,19X,'TF',/) 

DO 650 J=1,NF 

WRITE (2,640) FPOINT(J),TFPOINT(J) 

', 10X,F15.5) 

650 CONTINUE 

WRITE (2,660) 

660 FORMAT (’1',///,39X,'OUTPUT DATA',///,33X,'TIME',9X, 

C 'CONCENTRATION (CE)',/) 

DO 680 J=l,ITMAX 
T = J * TINCR 
WRITE (2,670) T,CE(J) 

670 FORMAT (25X,F15.5,' ’,10X,F15.5) 

680 CONTINUE 
RETURN 
END 



Appendix E 


Sample COMBALL Input File 
For Freundlich Parameters Determined by Excel 
and A-D Fractional Breakthrough Curve 


FREUNDLICH 

9.98927 

1.28623 

0.10581 

0.50000 

10.00000 

0.02500 

100 

100 

0.02500 

0.05000 

0.07500 

0.10000 

0.12500 

0.15000 

0.17500 

0.20000 

0.22500 

0.25000 

0.27500 

0.30000 

0.32500 

0.35000 

0.37500 

0.40000 

0.42500 

0.45000 

0.47500 

0.50000 

0.52500 

0.55000 

0.57500 

0.60000 

0.62500 

0.65000 

0.67500 

0.70000 

0.72500 

0.75000 

0.77500 

0.80000 

0.82500 

0.85000 

0.87500 

0.90000 

0.92500 

0.95000 

0.97500 

1.00000 

1.02500 

1.05000 

1.07500 

1.10000 

1.12500 

1.15000 

1.17500 

1.20000 

1.22500 

1.25000 

1.27500 

1.30000 

1.32500 

1.35000 

1.37500 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0.00004 

0.00013 

0.00040 

0.00107 

0.00258 

0.00561 

0.01119 

0.02061 

0.03536 

0.05692 

0.08654 

0.12498 

0.17235 

0.22803 

0.29068 

0.35840 

0.42896 

0.50000 

0.56931 

0.63497 

0.69550 

0.74991 

0.79767 

0.83869 

0.87318 

0.90165 

0.92471 

0.94308 

0.95748 

0.96859 

0.97706 

0.98342 

0.98813 


1.40000 

1.42500 

1.45000 
1.47500 
1.50000 
1.52500 
1.55000 
1.57500 
1.60000 
1.62500 
1.65000 
1.67500 
1.70000 
1.72500 
1.75000 
1.77500 
1.80000 
1.82500 
1.85000 
1.87500 
1.90000 
1.92500 
1.95000 
1.97500 
2.00000 
2.02500 
2.05000 
2.07500 
2.10000 
2.12500 
2.15000 
2.17500 
2.20000 
2.22500 
2.25000 
2.27500 
2.30000 
2.32500 
2.35000 
2.37500 

2.40000 

2.42500 

2.45000 
2.47500 

10.00000 


0.99159 

0.99409 

0.99589 

0.99716 

0.99805 

0.99868 

0.99911 

0.99940 

0.99960 

0.99974 

0.99983 

0.99989 

0.99993 

0.99995 

0.99997 

0.99998 

0.99999 

0.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 


COMBALL Output 

(Mass Balance Results Portion is Not Included) 
FREUNDLICH MULTIWELL INJECTION SIMULATION 

INPUT DATA 

FREUNDLICH PARAMETERS 


FREUNDLICH COEFFICIENT 

(B) = 

9.989 

FREUNDLICH EXPONENT (A) 

= 

1.286 

INJECTED CONCENTRATION 

(CO) = 

.106 

INJECTION TIME (DELT) 

= 

.500 


ISOTHERM 

A 

C = (B) (C) 

SORBED 


FRACTIONAL BREAKTHROUGH CURVE 


F 


TF 


.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00001 

.00004 

.00013 

.00040 

.00107 

.00258 

.00561 

.01119 

.02061 

.03536 

.05692 


.02500 

.05000 

.07500 

.10000 

.12500 

.15000 

.17500 

.20000 

.22500 

.25000 

.27500 

.30000 

.32500 

.35000 

.37500 

.40000 

.42500 

.45000 

.47500 

.50000 

.52500 

.55000 

.57500 

.60000 

.62500 

.65000 

.67500 

.70000 

.72500 

.75000 

.77500 

.80000 



.08654 

.12498 

.17235 

.22803 

.29068 

.35840 

.42896 

.50000 

1 

.56931 

1 

.63497 

1 

.69550 

1 

.74991 

1 

.79767 

1 

.83869 

1 

.87318 

1 

.90165 

1 

.92471 

1 

.94308 

1 

.95748 

1 

.96859 

1 

.97706 

1 

.98342 

1 

.98813 

1 

.99159 

1 

.99409 

1 

.99589 

1 

.99716 

1 

.99805 

1 

.99868 

1 

.99911 

1 

.99940 

1 

.99960 

1 

.99974 

1 

.99983 

1 

.99989 

1 

.99993 

1 

.99995 

1 

.99997 

1 

.99998 

1 

.99999 

1 

.99999 

1 

1.00000 

1 

1.00000 

1 . 

1.00000 

1 

1.00000 

1 , 

1.00000 

1 

1.00000 

1 , 

1.00000 

2 

1.00000 

2 

1.00000 

2 

1.00000 

2 

1.00000 

2 

1.00000 

2, 

1.00000 

2, 

1.00000 

2 

1.00000 

2, 

1.00000 

2, 

1.00000 

2 

1.00000 

2, 

1.00000 

2, 

1.00000 

2 . 

1.00000 

2. 

1.00000 

2 . 

1.00000 

2. 

1.00000 

2, 

1.00000 

2, 

1.00000 

2, 

1.00000 

10, 


.82500 

.85000 

.87500 

.90000 

.92500 

.95000 

.97500 

.00000 

.02500 

.05000 

.07500 

.10000 

.12500 

.15000 

.17500 

.20000 

.22500 

.25000 

.27500 

.30000 

.32500 

.35000 

.37500 

.40000 

.42500 

.45000 

.47500 

.50000 

.52500 

.55000 

.57500 

.60000 

.62500 

65000 

67500 

70000 

72500 

75000 

77500 

80000 

82500 

85000 

87500 

90000 

92500 

95000 

97500 

00000 

02500 

05000 

07500 

10000 

12500 

15000 

17500 

20000 

22500 

25000 

27500 

30000 

32500 

35000 

37500 

40000 

42500 

45000 

47500 

00000 


OUTPUT DATA 


TIME 


CONCENTRATION (CE) 


.02500 

.05000 

.07500 

.10000 

.12500 

.15000 

.17500 

.20000 

.22500 

.25000 

.27500 

.30000 

.32500 

.35000 

.37500 

.40000 

.42500 

.45000 

.47500 

.50000 

.52500 

.55000 

.57500 

.60000 

.62500 

.65000 

.67500 

.70000 

.72500 

.75000 

.77500 

.80000 

.82500 

.85000 

.87500 

.90000 

.92500 

.95000 

.97500 

1.00000 

1.02500 

1.05000 

1.07500 

1.10000 

1.12500 

1.15000 

1.17500 

1.20000 

1.22500 

1.25000 

1.27500 

1.30000 

1.32500 

1.35000 

1.37500 

1.40000 

1.42500 

1.45000 

1.47500 

1.50000 

1.52500 

1.55000 

1.57500 

1.60000 

1.62500 

1.65000 

1.67500 

1.70000 


.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00001 

.00001 

.00002 

.00003 

.00005 

.00007 

.00009 

.00012 

.00015 

.00019 

.00023 

.00026 

.00030 

.00034 

.00037 

.00040 

.00042 

.00045 

.00046 

.00048 

.00049 

.00050 

.00051 

.00052 

.00052 

.00053 

.00053 

.00054 

.00054 

.00054 

.00055 

.00055 

.00055 

.00056 

.00056 

.00057 

.00057 

.00058 

.00058 

.00059 



1.72500 

1.75000 

1.77500 
1.80000 

1.82500 

1.85000 

1.87500 

1.90000 

1.92500 

1.95000 

1.97500 
2.00000 
2.02500 
2.05000 
2.07500 
2.10000 

2.12500 

2.15000 

2.17500 

2.20000 

2.22500 

2.25000 

2.27500 

2.30000 

2.32500 

2.35000 

2.37500 

2.40000 

2.42500 

2.45000 

2.47500 

2.50000 
2.52500 
2.55000 
2.57500 
2.60000 
2.62500 
2.65000 
2.67500 
2.70000 

2.72500 

2.75000 

2.77500 
2.80000 

2.82500 

2.85000 

2.87500 

2.90000 

2.92500 

2.95000 

2.97500 
3.00000 
3.02500 
3.05000 
3.07500 
3.10000 

3.12500 

3.15000 

3.17500 
3.20000 

3.22500 

3.25000 

3.27500 

3.30000 

3.32500 

3.35000 

3.37500 

3.40000 

3.42500 

3.45000 

3.47500 

3.50000 


.00060 

.00061 

.00061 

.00062 

.00063 

.00064 

.00065 

.00066 

.00068 

.00069 

.00070 

.00072 

.00073 

.00075 

.00077 

.00079 

.00080 

.00083 

.00085 

.00087 

.00089 

.00092 

.00094 

.00097 

.00100 

.00103 

.00106 

.00109 

.00113 

.00116 

.00120 

.00124 

.00128 

.00132 

.00136 

.00140 

.00145 

.00149 

.00154 

.00158 

.00163 

.00168 

.00174 

.00179 

.00185 

.00191 

.00197 

.00204 

.00210 

.00217 

.00224 

.00232 

.00239 

.00247 

.00255 

.00263 

.00272 

.00281 

.00290 

.00299 

.00309 

.00319 

.00329 

.00340 

.00351 

.00362 

.00373 

.00385 

.00397 

.00409 

.00422 

.00435 



3.52500 

3.55000 

3.57500 

3.60000 

3.62500 

3.65000 

3.67500 

3.70000 

3.72500 

3.75000 

3.77500 

3.80000 

3.82500 

3.85000 

3.87500 

3.90000 

3.92500 

3.95000 

3.97500 
4.00000 
4.02500 
4.05000 
4.07500 
4.10000 

4.12500 

4.15000 

4.17500 

4.20000 

4.22500 

4.25000 

4.27500 

4.30000 
4.32500 
4.35000 
4.37500 
4.40000 
4.42500 
4.45000 
4.47500 
4.50000 

4.52500 

4.55000 

4.57500 

4.60000 

4.62500 

4.65000 

4.67500 

4.70000 

4.72500 

4.75000 

4.77500 

4.80000 

4.82500 

4.85000 

4.87500 

4.90000 

4.92500 

4.95000 

4.97500 
5.00000 
5.02500 
5.05000 
5.07500 
5.10000 

5.12500 

5.15000 

5.17500 
5.20000 

5.22500 

5.25000 

5.27500 

5.30000 


.00448 

.00462 

.00476 

.00491 

.00505 

.00520 

.00536 

.00552 

.00568 

.00584 

.00601 

.00619 

.00636 

.00654 

.00673 

.00692 

.00711 

.00731 

.00751 

.00771 

.00792 

.00814 

.00835 

.00858 

.00880 

.00904 

.00927 

.00951 

.00975 

.01000 

.01026 

.01051 

.01077 

.01104 

.01131 

.01158 

.01186 

.01214 

.01242 

.01271 

.01300 

.01330 

.01360 

.01389 

.01419 

.01450 

.01480 

.01511 

.01541 

.01571 

.01601 

.01631 

.01662 

.01693 

.01721 

.01749 

.01778 

.01807 

.01837 

.01862 

.01887 

.01912 

.01938 

.01964 

.01986 

.02006 

.02026 

.02046 

.02067 

.02086 

.02099 

.02112 


5.32500 

5.35000 

5.37500 

5.40000 

5.42500 

5.45000 

5.47500 

5.50000 

5.52500 

5.55000 

5.57500 

5.60000 

5.62500 

5.65000 

5.67500 

5.70000 

5.72500 

5.75000 

5.77500 
5 .'80000 

5.82500 

5.85000 

5.87500 

5.90000 

5.92500 

5.95000 

5.97500 
6.00000 
6.02500 
6.05000 
6.07500 
6.10000 
6.12500 
6.15000 
6.17500 
6.20000 
6.22500 
6.25000 
6.27500 
6.30000 

6.32500 

6.35000 

6.37500 

6.40000 

6.42500 

6.45000 

6.47500 

6.50000 

6.52500 

6.55000 

6.57500 

6.60000 

6.62500 

6.65000 

6.67500 

6.70000 

6.72500 

6.75000 

6.77500 
6.80000 

6.82500 

6.85000 

6.87500 

6.90000 

6.92500 

6.95000 

6.97500 
7.00000 
7.02500 
7.05000 
7.07500 
7.10000 


.02125 

.02138 

.02152 

.02158 

.02163 

.02167 

.02172 

.02177 

.02177 

.02173 

.02169 

.02164 

.02160 

.02153 

.02140 

.02127 

.02114 

.02101 

.02086 

.02066 

.02045 

.02024 

.02002 

.01980 

.01953 

.01926 

.01898 

.01869 

.01841 

.01809 

.01777 

.01744 

.01710 

.01677 

.01642 

.01606 

.01570 

.01534 

.01497 

.01460 

.01423 

.01386 

.01349 

.01310 

.01273 

.01237 

.01200 

.01163 

.01125 

.01088 

.01054 

.01019 

.00984 

.00948 

.00913 

.00882 

.00850 

.00817 

.00785 

.00753 

.00725 

.00696 

.00667 

.00638 

.00610 

.00585 

.00560 

.00535 

.00510 

.00485 

.00465 

.00444 



7.12500 

7.15000 

7.17500 

7.20000 

7.22500 

7.25000 

7.27500 

7.30000 

7.32500 

7.35000 

7.37500 

7.40000 

7.42500 

7.45000 

7.47500 

7.50000 

7.52500 

7.55000 

7.57500 

7.60000 

7.62500 

7.65000 

7.67500 

7.70000 

7.72500 

7.75000 

7.77500 

7.80000 

7.82500 

7.85000 

7.87500 

7.90000 
7.92500 
7.95000 
7.97500 
8.00000 
8.02500 
8.05000 
8.07500 
8.10000 

8.12500 

8.15000 

8.17500 

8.20000 

8.22500 

8.25000 

8.27500 

8.30000 

8.32500 

8.35000 

8.37500 

8.40000 

8.42500 

8.45000 

8.47500 

8.50000 

8.52500 

8.55000 

8.57500 

8.60000 

8.62500 

8.65000 

8.67500 

8.70000 

8.72500 

8.75000 

8.77500 

8.80000 

8.82500 

8.85000 

8.87500 

8.90000 


.00422 

.00401 

.00380 

.00363 

.00346 

.00328 

.00310 

.00293 

.00279 

.00265 

.00251 

.00236 

.00223 

.00212 

.00200 

.00189 

.00177 

.00167 

.00158 

.00149 

.00140 

.00131 

.00123 

.00116 

.00110 

.00103 

.00096 

.00090 

.00085 

.00079 

.00074 

.00069 

.00065 

.00061 

.00057 

.00053 

.00049 

.00046 

.00043 

.00040 

.00037 

.00034 

.00032 

.00030 

.00028 

.00026 

.00024 

.00022 

.00021 

.00019 

.00018 

.00016 

.00015 

.00014 

.00013 

.00012 

.00011 

.00010 

.00010 

.00009 

.00008 

.00008 

.00007 

.00006 

.00006 

.00005 

.52500 

1.22500 

1.47500 

1.70000 

1.92500 

2.10000 



8.92500 

8.95000 

8.97500 
9.00000 
9.02500 
9.05000 
9.07500 
9.10000 
9.12500 
9.15000 
9.17500 
9.20000 
9.22500 
9.25000 
9.27500 
9.30000 
9.32500 
9.35000 
9.37500 
9.40000 
9.42500 
9.45000 
9.47500 
9.50000 
9.52500 
9.55000 
9.57500 
9.60000 
9.62500 
9.65000 
9.67500 
9.70000 
9.72500 
9.75000 
9.77500 
9.80000 
9.82500 
9.85000 
9.87500 
9.90000 

9.92500 

9.95000 

9.97500 
10.00000 


2.25000 
2.37500 

2.50000 
2.60000 

2.67500 

2.75000 

2.82500 

2.90000 

2.95000 
3.00000 
3.07500 
3.12500 
3.15000 
3.20000 

3.25000 
3.27500 
3.32500 
3.35000 
3.40000 
3.42500 
3.45000 

3.50000 
3.52500 
3.55000 
3.57500 
3.60000 
3.62500 
3.65000 

3.67500 
3.70000 
3.72500 

3.75000 
3.77500 
3.80000 

3.82500 
3.85000 
3.87500 

3.90000 
3.92500 

3.95000 
3.97500 
4.00000 
4.02500 
4.05000 
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